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In this astro-particle lecture course I shall try to emphasize evidence of the new 
physics which we have in cosmological and astrophysical data. This includes 
support of the inflationary model, necessity of dark energy and of non-baryonic 
dark matter, the Grizen-Zatsepin-Kuzmin puzzle of the ultra-high energy cos- 
mic rays. 

1, INTRODUCTION 

The purpose of these CERN school lectures is to review the evidence for the new physics in cosmolog- 
ical and astrophysical data, and to give the minimal theoretical frameworks needed to understand and 
appreciate the evidence. Beyond any reasonable doubt, we have solid evidence for the new physics in 
these data. The strongest is the case for non-baryonic dark matter, followed by the case for dark energy. 
The possibility (though very speculative, since a consistent and working model has not been constructed 
yet) that the law of gravity should be changed instead is not excluded, but that would mean a new physics 
anyway. I will not engage in discussion of relevant particle physics model building, instead the reader 
is refereed to lectures by G. Gabadadze at this school |1 1. Another solid evidence for the new physics 
beyond the standard model is given by neutrino oscillations. I will not discuss this topic, it is covered in 
lectures by S. Petcov at this school |2|. The physics of cosmic rays is partially covered in lectures by A. 
Chilingaryan O, therefore, I restrict myself to the highest-energy part of the spectrum, which is related 
to Grizen-Zatsepin-Kuzmin puzzle and, possibly, to a new physics. 

There are many excellent reviews on the subject of Cosmology and Astroparticle physics, includ- 
ing lectures at previous CERN schools, for a recent one see |4, I've tried to be complimentary to 
these lectures to the extent it is possible, so many additional details can be found there. Proceedings of 
these schools can be found at http://physicschool.web.cern.ch/PhysicSchool Due to space and time limi- 
tations, I omit several very important traditional topics, most notably Big Bang Nucleosynthesis (see e.g. 
the review |7 1) and Baryogenesis (see e.g. the review 1 8 1). In the area covered, I've updated experimental 
results and resulting constraints, where applicable. The most important developments since the time of 
the previous school were: release of the first year observations of Cosmic Microwave Background Radia- 
tion (CMBR) by the Wilkinson Microwave Anisotropy Probe (WMAP) |9|, first data release by the Sloan 
Digital Sky Survey (SDSS) of three-dimensional distribution of galaxies 1 10 1, and the releasel 11 J of a 
statistically significant dataset of Supernovae la at large cosmological redshifts, z > I, which provide 
the first conclusive evidence for cosmic deceleration that preceded the current epoch of cosmic accel- 
eration. These are long awaited cosmological data of unprecedented quality, and with their appearance 
cosmology has truly entered the golden era and became a precision science. 

The plan of the lectures is as follows. In Section |3 I review the basics of cosmology: Friedman 
equations, Hubble expansion, cosmography. In Section|3 the Cosmic Microwave Background Radiation 
(CMBR) is discussed. In Section |^ I briefly review resent results on another cosmological probe - 
the large-scale distribution of galaxies. In Section |3 the evidence for the existence of dark energy is 
presented. S ections l6l and lOl review the evidence for a dark matter, and particle physics models of non- 
baryonic matter are briefly considered. In Section I review the basics of inflationary cosmology and 
discuss support of the inflationary model by the CMBR data. In Section|^ the physics of the Ultra-high 
energy cosmic rays is reviewed. 



2, BASICS OF COSMOLOGY 



2,1 Note on units and scales 

Length. Astronomers are measuring distances in parsecs, which is about 3.1 x 10^^ m or about 3.26 
light years, and is comparable with the distance to the closest star. PARSEC is an abbreviation for the 
distance to a star with a semi-annual PARallax of 1 arc SECond. The distance from our Sun to the 
Galactic center is 8 kiloparsecs, so the kpc is an appropriate unit when discussing galactic structure. The 
appropriate unit of extragalactic distance, however, is the megaparsec, or Mpc. The nearest large cluster 
of galaxies, the Virgo cluster, is at a distance of 20 Mpc. The size of the visible Universe is 4200 Mpc or 
13.7 billions of Ught years. 



Energy. Usually astronomers are measuring energy in ergs. E.g. the luminosity of our Sun is 4 x 
Iq33 gj.g while luminosity of bright quasars reaches 10^^ erg s~^. Galaxy like our Milky Way 
contains 10^^ stars, and there are 10^^ galaxies in the visible part of the Universe. Particle physicists 
are measuring energy in electron-volts, 1 erg = 6.2 x 10^^ eV, and usually are choosing units where the 
velocity of light and the Plank constant are set to unity, c = 1, h = 1, which I am using too, when 
convenient. In these units, for example, 1 Mpc = 1.6 x 10^^ eV~^. 



2.2 Dynamical Frameworks 

Dynamics is provided by General Relativity - the Einstein field equations 

Rf^u - \g,..R = SttG t^, , (1) 

where T^^ is a stress energy tensor describing the distribution of mass in space, G is Newton's gravi- 
tational constant, and the curvature R^^ is a complicated function of the metric and its first and second 
derivatives. Clearly, finding a general solution to a set of equations as complex as the Einstein field equa- 
tions is a hopeless task. The problem is simplified greatly considering mass distributions with special 
symmetries. The basic assumption underlying the construction of cosmological models is that of spatial 
homogeneity and isotropy. The most general space-time metric consistent with these symmetries is the 
Robertson- Walker metric: 

ds'^ = dt'^-a'^{t)dl'^ , (2) 

where a{t) is the dimensionless scale factor by which all distances vary as a function of cosmic time t. 
The scale factor contains all the dynamics of the Universe, while the vector product dl^ describes the 
geometry of the space. 



which can be either EucUdian, or positively or negatively curved. For the spatial 3-dimensional curvature 
we find, explicitly 

(\h ( k = -l Open 
(3)i?=4^ } k = Flat (3) 
^ A; = +1 Closed 



E.g., the space with A; = -1-1 can be thought of as a 3-dimensional sphere with a curvature being inversely 

proportional to the square of its radius. In this Section we will model the matter content of the Universe 
as a perfect fluid with energy density p and pressure p, for which the stress-energy tensor in the rest 
frame of the fluid is 
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With these assumptions the Einstein equations simpUfy to the Friedmann equations, which form the 
dynamical basis of cosmology 
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Let us have a look at the basic physics behind of these equations. 
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1. Differentiating Eq. ^ and subtracting Eq. Q we obtain 

f +3^(p + p) = 0, (7) 
at a 

which is nothing but energy-momentum conservation, 

t;., = o. (8) 

On the other hand, the result is nothing but the First Law of thermodynamics 

dE + pdV = TdS , (9) 

with dS = 0. Here E = pV = pa^Vo is energy, T is temperature and S is entropy of some (fixed) 
comoving volume Vq. Therefore, Friedmann expansion driven by an ideal fluid is isentropic, dS = 0. 
This is not unexpected, and relaxing the assumption of a perfect fluid will lead to entropy production. 
However, the dissipation is negligible in cosmological frameworks (except of special moments, like the 
initial matter creation and possible phase transitions, which will be considered separately) and isentropic 
expansion is a very good approximation. This gives a useful integral of the motion, S = const. On 
dimensional grounds, S oc T^a^Vb and we obtain frequently used relation between the scale factor and 
temperature in an expanding Universe 

a oc ^ . (10) 

To be precise, 

S = —g^T^a^ = const , (11) 
where the factor g^^ counts the effective number of relativistic degrees of freedom 

7 7 
9*= Yl + Yl 3j = {9B + ■^9f) ■ (12) 

j=bosons j=fermions 

At any given temperature, only particles with m <^ T should be counted, i. e. g^ is a function 
of temperature, which is shown in Fig. [2 For a gas of photons, g^ = 2. Considering the current epoch 
of the Universe evolution, we have to add the neutrino contribution, which will be discussed later and 
leads to a different account of effective degrees of freedom at temperatures below e+e^ annihilation in 
entropy, gs, and in energy density, gp. Namely, gs{To) = 3.909 and QpiTo) = 3.363. At temperatures 
above the electroweak scale g^ ~ 100 in the Standard Model. 

Let us give here also other thermodynamical relations, similar to Eq. (fTTl . but for the energy 
density, p, and particle number density, n 
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TTTT 9*T^ , (13) 
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n = ^{9b + -,9f)T\ (14) 
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Fig. 1: Number of relativistic degrees of freedom as a function of temperature. From Ref. 1121 . 



where C(3) = 1.202. These relations are a simple consequence of the integration of Bose-Einstein or 
Fermi-Dirac distributions 
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where q is particle momentum, the plus (minus) sign corresponds to fermions (bosons), and a = 1 in 
calculation of p, while a = in calculation of n (in the latter case the integral cannot be evaluated in 
terms of elementary functions and the Riemann (^-function appears). With the use of Eq. the entropy 
density, Eq. (fTTl . can be found as s = Ap/'iT, since for the relativistic particles p = p/?, regardless of 
spin. 

2. Friedmann equation, Eq. can be interpreted within Newtonian mechanics. Indeed, let us first 
re-arrange it as 

-a pa = . (16) 

2 3 2 ^ ' 

Now, it is easy to see that for r = aro, the Friedmann equation takes the form of energy conservation 

for test particles bounded in the gravitational potential created by mass M = ^ pr^ , 

1 , GM krl 

-r^ = ^. (17) 

2 r 2 ^ ' 

We see that the constant k, which determines the sign of spacial curvature in the language of General 
relativity, also determines the sign of the binding energy 

fc = +1 Binding energy is negative, 

the Universe will recollapse 
k = —\ Binding energy is positive, 

the Universe will expand forever 



Therefore, the case of zero spatial curvature, or zero binding energy. A; = 0, is special and corresponds 
to fine tuning between initial kinetic and potential energies. Setting A; = in Eq. this fine tuning can 




Fig. 2: Hubble Diagram. From Ref. 1141 . 



be expressed as p = pc, where the critical density is defined as 

The critical density is proportional to the square of another fundamental parameter, 

a 

H = -. (19) 

a 

The present value of this parameter is called the Hubble constant. It describes the rate of the Universe 
expansion, and can be related to observations in the following way. Consider two points with a fixed 
comoving distance tq between them (this means that points do not feel any other forces and do not par- 
ticipate in any other motion beyond general expansion of the Universe). The physical distance between 
points increases as r{t) = a(t)ro, and we can find the relative velocity as 

a 

V = cirQ = — aro = Hr . (20) 
a 

The relation v = Hr is called the Hubble law. This is shown in Fig. |2l The left panel is original 
data used by E. Hubble, the right panel presents recent data from Ref. 1 13|. We will discuss it in more 
detail later on, especially in relation to observations. But the first thing we may notice is that according 
to Hubble law, u ~ 1 at r ~ H^^. Separations (or wavelengths) of this order are therefore special 
in cosmology and mean super-horizon length scale. At smaller separations, Newtonian gravity should 
be valid. Einstein equations tell us that energy conservation in Newtonian mechanics, Eq. and 
the first law of thermodynamics, Eq. Q, applied to the Universe as a whole, can be extended beyond 
horizons without any change. The second Friedmann equation, Eq. can be derived as a consequence 
of these two equations. However, hardly would Newton do it, even if he had known the first law of 
thermodynamics. Indeed, in Eq. ^ we recognize the Newton's second law, F = mv with F being the 
gravity force, and energy conservation is derived from the equations of motion, not vice versa. 

We see that, according to Einstein's theory, the force law is modified. Not only mass gravitates, 
but the pressure, too, gives its contribution into gravitational force. This is a very important modification, 
since pressure can be negative, leading to anti-gravity and to accelerated expansion. As we will see, this 
stage of expansion may have lead to creation of the Universe in our classical understanding. At present. 



the Universe expansion seems to be dominated by anti-gravity as well. This has an interesting history. 
Newton did not know that one should worry about horizons, but he worried why the Universe does not 
collapse under the pull of gravity. Einstein was worried about this too. He added (1917) a cosmological 
constant to the equations of motion, thinking that it will make the Universe static. (As we will see, 
the cosmological constant corresponds to a vacuum with non-zero energy density and negative pressure, 
p = —p.) However, Friedmann had shown (1922) that the Universe will not be static anyway. After 
some debate, Einstein admitted his mistake and called the introduction of a cosmological constant "the 
greatest blunder of my life". 

So, why did the Universe not collapse under the pull of gravity? Resolution is in awkward initial 
conditions called the Big Bang, where velocity in Eq. dTTl is highly tuned to potential energy, leading 
to practically zero spatial curvature and to p = pc- This implies enormous fine-tuning for the Universe 
to survive till present. Such fine-tuning is hard to accept, and a modification of classical cosmology 
was called for. We will see how modern inflationary cosmology solves the problem of initial conditions. 
Again, the resolution is in anti-gravity caused by negative pressure. 

2.3 Matter content in the Universe 

To solve Friedmann equations, Eq. (Isll-©, one has to specify the matter content of the Universe and the 
equation of state for each of the constituents. To fit current observations we need at least four components 

- Radiation (relativistic degrees of freedom). Today this component consists of the photons and 
neutrino and gives negligible contribution into total energy density. However, it was major fraction 
at early times. 

- Baryonic matter. Makes up the observable world today. 

- Dark matter. Was not directly detected yet, but should be there. Constitutes major matter fraction 
today. Has rather long observational history and can be fitted within frameworks of modem particle 
physics nicely, at the price of "moderate" tuning of parameters to provide required fraction of 
matter. 

- Dark energy. Looks like it also should be there. It provides the major fraction of the total energy 
density today. Was not anticipated and appears as a biggest surprise and challenge for particle 
physics, though conceptually it can be very simple, being just a "cosmological constant" or vacuum 
energy. 

Equations of state. Each of these components have very simple equation of state, parameterized by a 
single constant w 

w = -. (21) 
P 

With a constant w, solutions of the first law of thermodynamics, Eq. Q, are readily found 

p{t) = a(t)-3(i+-) po , (22) 

where po stands for the present day density. For example, for ordinary forms of matter we have 

- Radiation: w = ^ and p = po- Result can be understood as a simple consequence of entropy 
conservation, aT = const, since for radiation p ocT^. 

- Matter: w = and p = a^^ pp. Result can be understood as a simple consequence of particle 
number conservation, na^ = const and p = mn, where m is particle mass. 

For hypothetical matter, which may play the role of dark energy, w is negative 

- Network of cosmic strings: w = and p = po. 

- Network of domain walls: w = and p = a^^ pq. 



- Cosmological constant: w = —1 and p = pQ. Result can be understood as a consequence of the 
Lorentz invariance of a vacuum, which restricts to be proportional to the Kronecker tensor, 
T// = V 5^. Comparing with Eq. @ we find p = —p, orw = —l. 

Law of expansion. Friedmann equation (|5ll in a spatially flat Universe and with a single matter com- 
ponent, which energy density evolves according to Eq. (l22l . has the solution 

a = (t/to)3{^ . (23) 

In particular 

- Radiation: w = a 

- Matter: w = 0, a 

- Cosmological constant: w 

2.4 Cosmological Parameters. 

These are used to parameterize the Fiedmann equation and its solution a{t). Let me first summarize the 
current knowledge of numerical values of those parameters which were introduced already; later in the 
course we will discuss how these values were deduced. 



Cosmological parameters 



Symbol & Definition 


Description 


Present value 


t 

H = a/a 

p^ = m'^/s^G 

Vt = p/pc 

^CMB — PcMbI Pc 
= Pb/Pc 
— Pm/ Pc 
= pa/ Pc 


Age of the Universe 
Hubble constant 
Critical density 
Omega 

Fraction of CMB photons 
Baryonic fraction 
Matter fraction 
Dark Energy fraction 


to = (13.7 ± 0.2) Gyr 
Hq = 71 km s"^ Mpc~^ 
Pc = 10 h/^ GeV 

= 1.02 ±0.02 
J^cMB = 2.4 • 10-5/^-2 

= 0.044 ± 0.004 
Q.m = 0.27 ±0.04 
VLa = 0.73 ± 0.04 



2.5 Cosmography. 

We can define cosmography (this is my custom definition for these lectures) as a part of cosmology which 
tries to map observations into reconstruction of the scale factor. I.e., the goal is to find and tabulate the 
function a{t). This is important in many respects: e.g. it allows to determine the matter content in 
the Universe (assuming the Friedmann equations are correct). One simple and straightforward way of 
tabulating the function is in determining its coefficients in Taylor expansion. This can be done making 
a Taylor decomposition around present time, t = t^. The value of the scale factor at any moment of 
time can be fixed arbitrarily, we can use this freedom to choose a (to) = 1- The second term in the 
Tailor decomposition is naturally the value of the Hubble constant, Hq, Eq. (fT9l . It gives us the rate 
of the Universe expansion at present and can be measured using the Hubble law, Eq. (l20t . One can go 
further in this decomposition and define the second derivative of the universe at present, the "deceleration 
parameter", and so on. We will not do it (at least at this point), since modern observations probe the whole 
function a{t). Therefore, let us start with the preparation of the necessary machinery which allows us to 
deduce a{t) from observations. 

Let me stress now that Eq. d20b involves some degree of cheating since it is not a relation between 
the observables. (However, it is a valid relation for small separations.) To apply the Hubble law to 
observations, we have to derive its generalization, which would connect quantities we can measure. 



= — 1. This case is special, and a = e 



Redshift. Looking at distant objects we see only the light they emit. How can physical quantities like 
distance and velocity be derived? Recall how a police officer determines the speed of a car. A similar 
principle is used in cosmology to determine the velocity of distant bodies. The shift of emission lines 
with respect to the frequency measurements by the local observer is related to velocity, and is used as 
an observable instead of the velocity. Systematic recession of objects, or cosmological expansion, leads 
to redshift. Note that cosmological redshift is not entirely due to the Doppler effect, but, rather, can be 
interpreted as a mixture of the Doppler effect and of the gravitational redshift. 

Let us relate now the redshift to cosmological expansion, aja. To this end, we consider photon 
trajectories in a cosmological background with metric Eq. The trajectory is given by ds^ = 0. 
Since the overall scale factor does not change the solutions of ds^ = 0, it is convenient to introduce the 
conformal time rj defined as 

dt = a{7]) dr] . (24) 

It is sufficient to consider radial trajectories with the observer at the center, and I restrict myself to a 
spatially flat metric ds'^ = o?{drf' — dx^) = 0, where x denotes the radial coordinate. The solution 
of di]'^ — d\^ = is X = + const. Since the comoving distance between source and observer does 
not change, the conformal time interval between two light pulses is the same at the point of emission and 
at the point of observation, Ar/ = const. Using the definition of conformal time, drj = dt/a, we find 

At _ 

emission — detection ■ 

a a 

Therefore, for a signal frequency we get cj^a^ = WeOe- Defining (measurable) redshift as 



we obtain 



(25) 



l + z=^. (26) 

It is convenient to normalize the scale factor by the condition = 1 at the point of detection, and to 
consider the scale factor at the point of emission as a function of redshift z. On the basis of this relation, 
the expansion history of the Universe can be parameterized as 

a(A = ^ . (27) 

The differential form of Eq. d27b is da/ dz = —cP'. For future use, let us find now the relation between 

di] and dz 

dri dt da , a , dz 

^^ = :^:r:r^^= -^^^= - tttt • ^^s) 

dt da dz a -"(-^j 

Observing that d-q = —dx, we obtain the Hubble low for small separations, dz = Hdx- At this point, 
we have succeeded in replacing the velocity by the redshift. Now we aim to relate the distance to some 
other quantity, directly measurable in cosmology. 



Luminosity distance. Looking at distant objects we see only the light they emit. How can physical 
quantities like distance and velocity be derived? There are several ways to introduce a quantity related 
to the distance: different definitions are not equivalent in curved space-time. A definition based on 
flux measurements is the appropriate one, if "standard candles" can be found and used. Detected flux 
[erg s^^ sm^^] is inversely proportional to the distance from a source, F oc D^^. Namely, if L is 
intrinsic luminosity [erg s^^], we have 

^2 L 



D^h called the "luminosity distance". For this technique to work, one has to find a set of sources with 
a known or calibrated luminosity. If such can be defined, they are called standard candles. 

To see how the luminosity distance enters the Hubble law, let us consider a space-time with the 
metric ds'^ = {drf — dy^ — \^ dQ). Now, go through the following list to find out what happens with 
the flux emitted into a frequency interval dv by a source located at redshif z: 

- Surface area at the point of detection is 47rx^. (Recall our choice = 1). 

- Energy and arrival rates are redshifted between the points of emission and detection. This reduces 
the flux by {l + zf. 

- Opposing this tendency, the bandwidth dv is reduced by (1 + z). 

- Photons observed at a frequency u, were emitted at [1 + z) v. 
Therefore, the measured spectral flux density is 

5(.) = Ml±il4 . (30) 

For the bolometric flux (i.e. integrated over z/) we find 

^=A 27T— (31) 

Comparing this with the definition, Eq. ( l29t . we find for the luminosity distance 

Dl = {1 + z)x, (32) 
where x is the comoving distance between the point of emission and the point of detection 

In the last equaUty we have used Eq. d28l . Therefore, the generahzation of the Hubble law, which can be 
used in observational cosmology, can be written as 



{l + z)x{z) = ^J^ . (34) 

Parameterization of H(z). Let us express now the function H{z) in the r.h.s. of Eq. through the 
cosmological parameters. First, we define the ratio of the total energy density to the critical one 

n = ^ . (35) 

Pc 

The present day value is referred to as Qq. Similarly, for each energy component we denote its present 
day fractional contribution as Qi = Pi/pc- With these definitions, the Friedmann equation Q for a 
spatially flat Universe can be re-written as 

i 

or 

= Hi ^ 17, (1 + 2)3(1+^0 . (36) 

i 

Here I have used 

Pi — Pifl a — Pc i'i U + 5 

and expressed the scale factor as a function of redshif t according to Eq. (126 1 . and used the definition of 
the critical density, = SvrG/Oc/S. 

Parameterization d36l is ready for use in Eq. (l33l for the comoving distance. In particular, this 
finalizes expression Eq. (l34l for the luminosity distance. 



3. Cosmic Microwave Background Radiation 

The Universe is filled with radiation which is left-over from the Big Bang. The name for this first light 
is Cosmic Microwave Background Radiation (CMBR). Measurements of tiny fluctuations (anisotropy) 
in CMBR temperature give a wealth of cosmological information and became a most powerful probe of 
cosmology. 

This radiation was predicted by Georgi Gamov in 1946, who estimated its temperature to be 
~ 5 K°. Gamov was trying to understand the origin of chemical elements and their abundances. Most 
abundant, after hydrogen, is helium, with its shear being ~ 25%. One possibility which Gamov consid- 
ered was nucleo-synthesis of He out of H in stars. Dividing the total integrated luminosity of the stars by 
the energy released in one reaction, he estimated the number of produced He nuclei. This number was 
too small in comparison with observations. Gamov assumed then the oven where the light elements were 
cooked-up was the hot early Universe. He calculated abundances of elements successfully and found that 
the redshifted relic of thermal radiation left over from this hot early epoch should correspond to ~ 5 K° 
at present. In one stroke G. Gamov founded two pillars (out of four) on which modem cosmology rests: 
CMBR and Big Bang Nucleosynthesis (BBN). Hot Big Bang was bom. 

Cosmic microwave background was accidentally discovered by Penzias and Wilson lITSl at Bell 
Labs in 1965 as the excess antenna temperature which, within the limits of their observations, was 
isotropic, unpolarized, and free from seasonal variations. A possible explanation for the observed excess 
noise temperature was immediately given by Dicke, Peebles, Roll, and Wilkinson and was published in a 
companion letter in the same issue 1 16|. They were preparing dedicated search experiment, but were one 
month late. Penzias and Wilson measured the excess temperature as ~ 3.5 ± 1 K° . It is interesting to 
note that the first (unrecognized) direct measurements of the CMB radiation was done by T. Shmaonov 
at Pulkovo in 1955, also as an excess noise while calibrating the RATAN antenna |17|. He published 
the temperature as (3.7 it 3.7) K° . Prior to this, in 1940, Andrew McKellar |18| had observed the 
population of excited rotational states of CN molecules in interstellar absorption lines, concluding that it 
was consistent with being in thermal equilibrium with a temperature of 2.7 K°. Its significance was 
unappreciated and the result essentially forgotten. Finally, before the discovery, in 1964 Doroshkevich 
and Novikov in an unnoticed paper emphasized [191 the detectability of a microwave blackbody as a basic 
test of Gamov's Hot Big Bang model. To me, as to theorist, the detection of CMBR looks nowadays like 
an easy problem. Indeed, a few percent of the "snow" on TV screens is due to CMBR. 

The spectrum of CMBR is a perfect blackbody, with a temperature T = 2.725 ± 0.002 K° as 
measured by modern instruments. This corresponds to 410.4 photons per cubic centimeter or to the flux 
of 10 trillion photons per second per squared centimeter. The temperature is slightly different in different 
patches of the sky - to 1 part in 100,000. And this is most important: the spectrum of this tiny fluctuations 
tells us a lot about the fundamental properties of the Universe. 

CMBR is the oldest light in the Universe. When registering it, we are looking directly at the 
deepest past we can, using photons. These photons had traveled the longest distances without being 
affected by scattering, and geometrically came out almost from the universe Horizon. More precisely, 
the CMB comes from the surface of the last scattering. We cannot see past this surface. That is because 
at early times the Universe was ionized and not transparent for radiation. With expansion, it cooled 
down, and when hydrogen recombined, the universe became transparent. Therefore the CMBR gives us 
a snapshot of the baby Universe at this time, which is called the time of last scattering. Let us determine 
when the last scattering had occurred in the early Universe. 

3.1 Hydrogen recombination 

At temperatures greater than a few thousands K, the ionized plasma in the Universe consisted of mostly 
protons, electrons, and photons, with a small fraction of helium nuclei and a tiny trace of some other 
light elements. To a good approximation we can consider only the hydrogen. Matter is then ionized at 



temperatures higher than the hydrogen ionization energy E'ion = 13.6 eV. At lower T neutral atoms start 
to form. The baryonic matter is in thermal equilibrium and the equilibrium fraction of ionized hydrogen 
can be described by the Saha equation (see Ref. |4| for more details) 

UhLy'e-E-^oJT ^ (37) 



rij^ \ 2-ir 

here rig, Up and Ujj are the number densities of electrons, protons, and neutral hydrogen respectively. 
Plasma is electrically neutral, i.e. rig = Up. To find the closed relation for the fraction of ionized atoms, 
X = np/{np + n^) = Up/n^, we need the relation between the baryon number density, n^, and 
temperature. This relation can be parameterized with the help of an important cosmological parameter 
called baryon asymmetry 

ri='^ = !hi±^ = (6.1 ± 0.3) X 10-1° , (38) 
where is the number density of photons 

n, = , (39) 

and (^(3) = 1.202, see Eq. (IT4l . Baryon asymmetry can be estimated by an order of magnitude by simply 
counting the number of baryons, rj = 2.68 x IQ-^ilfc/i^. This is not the most precise method, though; 
the value presented in Eq. (l38l was obtained from fitting the spectrum of CMBR fluctuations, see below. 
Nowadays, this is the most precise baryometer. Prior to this, the best estimates were obtained comparing 
BEN predictions of element abundances to observations. Defining recombination as the temperature 
when X = 0.1, we find Tree « 0.3eV. 

The Universe became transparent for radiation when the mean free path of photons became com- 
parable to the size of the Universe at that time. Photons scatter mainly on electrons and we find that the 
Universe became transparent when 

{a^ene)-^ ~ t . (40) 

Here, cr^e = 87ra2/3mg is the Compton cross-section. For the temperature of last scattering we find 
Tis ~ 0.26 eV. Taking the ratio to the current CMBR temeperature we find z\s ~ 1000. 



3.2 Spectrum is not distorted by red-sliift 

Prior to recombination photons were in thermal equilibrium. Therefore, at last scattering they have 
Planck spectrum 

""^^^ = exp(ii;i,/ri,) - 1 ■ 

Since then, particle momenta are red-shifted, p = k/a. Since photons are massless, E = p, their energies 
are red-shifted at the same rate, EqUo = -BisOis, and the spectrum becomes 

1 1 

exp(£^o/ais7is) - 1 exp(^o/7o) - 1 ' 

where we have used the notation To = ais^is. Therefore, after decoupling, the shape of the spectrum is 
not distorted. This statement would not be true for massive particles, E'^ = {p/a)^ + m?. 

Measuring CMBR, we should still see the Planckian spectrum, but with red-shifted temperature. 
Clearly, this conclusion is true not only for cosmological red-shift, but for the gravitational red-shifts as 
well. E.g., fluctuations in the gravitational potential at the last scattering surface should cause fluctuations 
in CMBR temperature, but do not distort the spectrum. 



Fig. 3: Left panel: uniform spectrum; error-bars are a small fraction of the line thickness. Right panel: dipole spectrum; vertical 
lines indicate one a uncertainties. From Ref. 1201 . 



3.3 Dipole spectrum 

Intensity of the CMB radiation is a function of the frequency 27rz^ = E and the direction on the sky (/, 6). 
As a function of 6) it can be decomposed in spherical harmonics. Coefficients will be the functions of 
V. The first two terms in this decomposition are 



where D{l,b) = cos 6, and 9 is an angle between observation and the maximum of the dipole Iq = 
263.85°, bo = 48.25°. The dipole is caused by our motion with respect to CMBR (which is composed 
of the motion of the Sun in the Galaxy and the Galaxy's own motion in the Local Cluster of galaxies). 
It gives the direction of this motion, Iq, bo, which roughly coincides with the direction towards Virgo. 
The dipole induced by the velocity v is vT cos 6. This gives the magnitude of Sun's peculiar velocity, 
(371 ± 1) km s-^. 

Let X = E/T. The monopole term should have the usual black-body spectrum lo^iy) oc / (e^ — 1 ) . 
The dipole spectrum is actually distorted, because the Doppler frequency shift depends upon direction. 
The dipole spectrum can be found as a term linear in v in the Taylor decomposition of S{v, I, b), with 
the result d{v) oc /(e^ —1)^; for a recent discussion see Ref. |21|. Functions and d{v) are 
shown in Fig. |3j left and right panels respectively. Both agree with theoretical expectation. 

3.4 Multipoles 

Monopole and dipole contributions to CMBR, Eq. dTO . can be subtracted. Emission of our Galaxy 
and various extragalactic sources can be subtracted also. This procedure uses the fact that the relic 
CMBR signal has a black body spectrum , which allows to distinguish it from other forms of radiation: 
measurements of the intensity at different frequencies allow to subtract contaminating foregrounds. What 
remains corresponds to the primordial cosmological pattern of temperature fluctuations, which is shown 
in Fig.|4] The upper panel presents the results of early COBE experiments [22 1, the lower panel shows 
the results of a recent WMAP experiment [9 1. 

The temperature anisotropy, T(n), as a function of viewing direction vector n, is naturally ex- 
panded in a spherical harmonic basis, Yij^ 



S{u,l,b) = Io{u) + D{l,b) d{iy) + . . . 



(41) 
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Fig. 4: Pattern of primordial temperature fluctuations in Galactic coordinates, from Ref. f9\. The WMAP map has 30 times 
finer resolution than the COBE map. 



The coefficients in this decomposition, aim, define the angular power spectrum, Ci 
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(43) 
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The CMBR angular power spectrum as measured by WMAP experiment is shown in Fig. |5] The har- 
monic index / is related to the angular scale 6* as ? 7r/6', so the first peak, at I « 220, would correspond 
to an angular scale of about one degree. Assuming random phases, the r.m.s. temperature fluctuation 
assosiated with the angular scale / can be found as 



Another representation of temperature fluctuations is given by the angular correlation function, which is 
related to Ci as 



where Pi is the Legendre polynomial of order /. 
3.5 Tool of Precision Cosmology 

The functional form of the CMBR power spectrum is very sensitive to both the various cosmological 
parameters and to the shape, strength and nature of primordial fluctuations. Measurements of the power 
spectrum provide us with a wealth of cosmological information at an unprecedented level of precision. 

Right after the discovery of CMBR, it was realized that fluctuations in its temperature should 
have fundamental significance as a reflection of the seed perturbations which grew into galaxies and 
clusters. In a pure baryonic Universe it was expected that the level of fluctuations should be of the order 
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Fig. 5: The angular power spectrum of primordial CMBR temperature fluctuations is shown in the upper panel, from Ref. 
Black symbols (I < 700) are WMAP measurements, data points at smaller angular scales represent CBI and ACBAR experi- 
ments. Lower panel shows the temperature-polarization cross-power spectrum. 

ST/T ~ 10~^ — 10~^. Mesurements of the CMBR anisotropy with ever-increasing accuracy have begun. 
Once the temperature fluctuations were shown to be less than one part in a thousand, it became clear that 
baryonic density fluctuations did not have time to evolve into the nonlinear structures visible today. A 
gravitationally dominant dark matter component was invoked. Eventually, fluctuations were detected at 
the level of 6T/T ~ 10~^ [ 23 ], consistent with structure formation in Cold Dark Matter models with the 
Harrison-Zel'dovich spectrum of primordial perturbations motivated by cosmological inflation. Already 
this magnitude of 6T/T is very restrictive by itself. A partial set of best fit cosmological parameters, as 
derived from the recent measurments of CMBR anisotropics, is presented in Table I. 

The foundations of the theory of CMBR anisotropy were set out by Sachs & Wolfe l24l . Silk 
f25l, Peebles & Yu |26|, Syunyaev & Zel'Dovich |27|. The measured spectrum of CMBR power has a 
characteristic shape of multiple peaks. Positions of these peaks and their relative amplitudes aie sensitive 
to many cosmological parameters in a non-trivial way. Fitting the data to model predictions gives very 
accurate values for many of these parameters (though there are some degeneracies between deferent sets). 
Numerical calculations for different models were done already in Ref. f28 1, and power spectra exhibiting 
acoustic peaks (similar to those in Fig. |5ll were presented. It was realized, in particular, that positions of 
the peaks are shifted with respect to each other for adiabatic and isentropic primordial fluctuations. 

3.6 Acoustic oscillations 

Let us give a qualitative picture of why the CMBR power spectrum has a specific shape of a sequence of 
peaks, and explain how it depends on the values of particular cosmological parameters. Insight, sufficient 
for the purposes of these lectures, can be gained with the idealization of a perfect radiation fluid. In 
complete treatment, one has to follow the evolution of coupled radiation and metric fluctuations, i.e. to 
solve the linearized Einstein equations. However, essential physics of radiation (or matter) fluctuations 



can be extracted without going into the tedious algebra of General Relativity. It is sufficient to consider 
the energy-momentum conservation, Eq. To solve for metric perturbations, full treatment based on 
Einstein equations, Eq. Q, is needed of course. We will not do that here, but simply quote results for the 
evolution of the gravitational potentials (coincident in some important cases with the solutions for the 
Newtonian potentials). 

Perturbations of the ideal radiation fluid, p = p/3, can be separated into perturbations of its tem- 
perature, velocity and gravitational potential. In the general-relativistic treatment gravitational potential 
appears as a fractional perturbation of the scale factor in the perturbed metric 

ds'^ = a^{7]) [(1 + 2^)^772 _ (1 _ 2^)dx'dxj] . (46) 

Two equations contained in the energy-momentum conservation, T^'^ = (i-e. temporal /i = and 
spatial fi = i parts of this equation), written in metric d46l . can be combined to exclude the velocity 
perturbations. The resulting expression is simple 

ek + Y(^k = -Y^k + ^k. (47) 

Note that this equation is the exact result for a pure radiation fluid. Here, 9^ are Fourier amplitudes of 
6T/T with wavenumber k, and is a Fourier transform of gravitational potential. Analysis of solutions 
of the Einstein equations for ^ shows that = const in two important cases: 

1. For superhorizon scales, which are defined as krj <^ 1. 

2. For all scales in the case of matter dominated expansion, p = 0. 

In these situations the last term in the r.h.s. of Eq. Wfl . namely, can be neglected. The Einstein 
equations also restrict the initial conditions for fluctuations. For the adiabatic mode in the limit kr] <^ 1 
one finds 

Si = -2<P^, (48) 

where 6i = Sp/p. The adiabatic mode is defined as a perturbation in the total energy density. For the 
one component fluid, which we consider here, only the adiabatic mode can exist. Note that fractional 
perturbation of the scale factor in metric (146 1 . a(?7,x) = a{r]) + 5a(?y,x) = a{7]){l — <I>), can be 
expressed as perturbation of spatial curvature, see Eq. Therefore, adiabatic perturbations are also 
called curvature perturbations. Let us re-write Eq. (l48l for temperature perturbations: 

• Radiation domination, 5 = A5T /T, and we find 

= - Y • (49) 

• Matter domination, 6 = 36T/T, and we find 

Oi = . (50) 

Recall now that in the limit k-q <^ I the gravitational potentail is time-independent, $ = const. There- 
fore, Eq. (l47l has to be supplemented by the following initial conditions in the case of the adiabatic 
mode: 

9i^0, ei = 0. (51) 

Temperature fluctuations on largest scales. Let us consider the modes which had entered the horizon 
after matter-radiation equality, fcr/eq < 1- For those modes, <!> = all the way from initial moments till 
present, and the solutions of Eq. Wfl with adiabatic inital conditions is 

e{r]) + $ = (^. + $) COS (^-^^ . (52) 




Fig. 6: Sensitivity of the CMBR angular power spectrum to four fundamental cosmological parameters (a) the curvature as 
quantified by fitot (b) the dark energy as quantified by the cosmological constant Qa (wa = — 1) (c) the physical baryon 
density Qth^ (d) the physical matter density Qmh^, all varied around a fiducial model of fitot = 1, = 0.65, Qi,h^ = 0.02, 
fl^h^ = 0.147. FromRef |29|. 



As gravity tries to compress the fluid, radiation pressure resists resulting in acoustic oscillations. It is 
important that oscillations are synchronized. All modes have the same phase regardless of k. This is a 
consequence of 6i = 0, which is valid for all k. At the last scattering, the universe becomes transparent 
for the radiation and we see a snapshot of these oscillations at-q = t]-^^. 

To get its way to the observer, the radiation has to climb out of the gravitational wells, which 
are formed at the last scattering surface. Therefore the observed temperature fluctuations are 6'obs = 
9{r]J + or 

^obs = COS (^^^ , (53) 

where we have used Eq. dSOl . which relates initial values of 9 and ^. Note that overdense regions 
correspond to cold spots in the temperature map on the sky, since the gravitational potential is negative. 
This is the famous Sachs-Wolfe effect l24l . 



Acoustic peaks in CMBR. Modes caught in the extrema of their oscillation, A;n%/\/3 = nvr, will 
have enhanced fluctuations, yielding a fundamental scale, or frequency, related to the universe sound 
horizon, s* = rj^^/VS. By using a simple geometrical projection, this becomes an angular scale on the 
observed sky. In a spatially flat Universe, the position of the first peak corresponds to li 200, see 
below. Both minima and maxima of the cosine in Eq. d53b give peaks in the CMBR power spectrum, 
which follow a harmonic relationship, A:„ = nvr/s^, see Fig.|5l 

The amplitudes of the acoustic peaks are recovered correctly after the following effects are taken 
into account: 1) baryon loading; 2) time-dependence of <I> after horizon crossing in radiation dominated 



universe; 3) dissipation. 

The effect of baryons is exactly the same for the oscillator equation Eq. d52l . as if we had increased 
the mass of a load connected to a spring and oscillating in a constant gravitational field starting on top 
of an uncompressed coil at rest. The addition of baryons makes a deeper compressional phase, and 
therefore increases every other peak in the CMBR power spectrum. (First, third, fifth, . . . ) The CMBR 
power spectrum is a precise baryometer. 

Gravitational potentials are not constant, but decay inside the horizon during radiation domination. 
This decay drives the oscillations: it is timed to leave compressed fluid with no gravitational force to fight 
when the fluid turns around. Therefore, the amplitudes of the acoustic peaks increase as the cold dark 
matter fraction decreases, which allows to measure ilm- 

Dissipation leads to dumping of higher order peaks in the CMBR power spectrum. 

The dependence of the CMBR angular power spectrum on different cosmological parameters is 
shown in Fig. |6l 



Position of the first peak. Position of the first peak is determined by the angular size of the sound 
horizon at last scattering. Let us calculate here a similar quantity: the causal horizon (which is larger 
by a factor of \/3 in comparison with the sound horizon). The comoving distance traveled by light, 
ds^ = 0, from the "Big Bang" to redshift z is determined by a relation similar to Eq. (l33l . but with 
different integration limits 



r]{z) 



dz' 



(54) 



H{z') ' 

where H{z) is given by Eq. One has to integrate this relation with a complete set of Oj, but for 
simplicity let us consider here the Universe dominated by a single component pj 



viz) 



(1 + ^) 



7j ^0 



where jj = (1 + 3wj)/2. From the last scattering to z ~ 1, the Universe was matter dominated. 
Therefore, the causal horizon in matter dominated Universe {wj = 0) 

o 

ri{z) 



should give a reasonable first approximation. 

Consider two light rays registered at z = which were separated by a comoving distance x = vi^) 
at the moment of emission. Since both propagate in the metric ds^ = a^{d7f — dx^ — x^dO"^) = 0, we 
find for the angular size of horizon at last scattering 



Oh 



??(o) ^/rn; 




(55) 



• This gives the position of the first acoustic peak, / « 200. 

• Tells us that there were 10'* causally disconnected regions at the surface of last scattering. 



Horizon problem. Regions separated by more than > 2° have not been in the causal contact prior to 
the last scattering in the standard Friedmann cosmology. The microwave sky should not be homogeneous 
on scales > 2°. Yet, CMB is isotropic to better than 10~^ on all scales. Observations tell us that all 
modes were, indeed, synchronized according to adiabatic initial conditions, Eq. dSTT i. with only small 
initial perturbations present, <I>i <^ 1. This constitutes the so-called "Horizon problem" of standard 
cosmology. In Section0we will see how this problem is solved in frameworks of inflationary cosmology. 




Fig. 7: Left panel: comparison of CMB power spectra in the models with adiabatic and isocurvature initial perturbations, from 
Ref. I.30J. Right panel: adiabatic power spectra in comparison with spectra appearing in models seeded by topological defects, 
from Ref. 1311 . In this panel some older, pre-WMAP, data are also shown. 



Non-adiabatic perturbations. For the isocurvature perturbations, instead of Eq. dSTT i. the initial con- 
ditions are given by 

5i = 0, 5i / . (56) 

That is because, in this case, perturbation in total density (and therefore in curvature) are zero initially. 
As a consequence, in Eq. ( I52t we will have sine instead of cosine. Acoustic peaks will be shifted by half 
a period, see Fig.0 Therefore, isocurvature perturbations are ruled out by modern CMBR experiments. 

If density perturbations would be seeded by topological defects (e.g. cosmic strings), both sine and 
cosine will be present in the solution for temperature fluctuations, Eq. (l52t . That is because the source 
for is active inside the horizon and phases of 9^ will be random. Acoustic peaks will be absent, 
see Fig. Structure formation seeded primarily by topological defects is ruled out by modern CMBR 
experiments. 



4. Large scale distribution of galaxies 

Primordial cosmological fluctuations leave their imprint as CMBR anisotropics (discussed in the previ- 
ous Section), and as density perturbations which give rise to galaxies and clusters of galaxies. CMBR 
anisotropics are observed on a two dimensional surface of last scattering, and therefore are measured as 
a two dimensional power spectrum. On the other hand, the distribution of galaxies can be measured in 
three dimensions. (Two angular coordinates of the line of sight to a galaxy and its redshift.) Different 
physical processes influence the initial perturbations until they are transformed into CMBR fluctuations 
or fluctuations of the distribution of galaxies. This influence can be encoded as a function of momenta, 
the transfer function T{k), which simply maps the power spectrum of the initial perturbations into the 
observed power spectrum, and is a function of cosmological parameters. Therefore, the distribution of 
galaxies gives complimentary information with respect to CMBR anisotropics and helps to break degen- 
eracy between cosmological parameters and the initial spectrum. 

This is illustrated in Fig. [8l with CMBR data from WMAP and large-scale structure data from 
SDSS. The left panel corresponds to 95% constraints in the (ilm, h) plane. The shaded dark red region is 
ruled out by WMAP alone leaving the long banana region. This shows that these two basic cosmological 
parameters are not well constrained by WMAP alone. The shaded light red region is ruled out when 
SDSS information is added. The small (shown as white) region remain allowed. Note that the allowed 
region is in good agreement with a completely independent measurement by HST key project based on 





Fig. 8: Combined CMBR and large-scale structure constraints. Left panel: 95% constraints in the (r2m,/i) plane. Right 
panel: 95% constraints in the (fim, CIa) plane. From Ref. 1101 . 



entirely different physics. The combined WMAP + SDSS constraint is even tighter than HST project 
measurement. 

One should bear in mind that there are caveats here. In deriving the WMAP+SDSS constraints 
which are shown in this figure, it was assumed that the universe is spatially flat, neutrino have negligible 
masses and the primordial spectrum is a pure power law. Without these priors the constraints are less 
tight. 

The constraints in the {^Im, ^a) plane with the assumption about spatial curvature being relaxed 
is shown in Fig. [8l right panel. The shaded dark red region is ruled out by WMAP alone, illustrating 
the well-known geometric degeneracy between models that all have the same acoustic peak locations. 
The shaded light red region is ruled out when adding SDSS information. Continuing inwards, the next 
two regions are ruled out allowing the assumption that re-ionization optical depth r < 0.3 and when 
supernova SN la information is included. Qa > is required with high confidence only when CMBR is 
combined with galaxy clustering information, or SN la information, see the next Section. 

5. Dark energy 

Something which is often called "Dark Energy" reveals itself in a variety of cosmological and astro- 
physical observations. This form of matter gravitates, but does not cluster. Contrary to radiation or dark 
matter, the dark energy causes the accelerated expansion of the universe. The need for it was hinted 
long ago to resolve conflict between the measured Hubble constant and the lower limits on the age of 
the Universe. Without the cosmological constant it was also not possible to obtain the correct growth of 
large scale structures in the Q = 1 Universe. Recently, the presence of dark energy was derived from the 
spectrum CMBR anisotropics and directly detected in the Hubble diagram of high redshift supemovae. 

Age of the Universe. If there is no dark energy, the Universe should be matter dominated and should 
expand according to a = {t/to)'^/^. Differentiating this expansion law we find 

Ht = 2/3 . (57) 

The value of the Hubble constant, as derived by the Hubble Key Project from the Hubble diagram, is 
72 ± 8 km s-i Mpc-i |14|, see Fig. |2l On the other hand, the lower bound on the age of the Universe 
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can be established estimating the ages of various objects it consists of. For example, the temperature of 
the coldest white dwarfs in globular clusters yields a cluster age of 12.7 ± 0.7 Gyr (T3\ . This gives 

i?oio > 0.93 ±0.12 , (58) 

in clear disagreement with Eq. d58l . In other words, the Universe appears much younger than the ages 
of the oldest objects in it. The critical density Universe, Q = 1, cannot consist of pressureless matter if 
measurements of the Hubble constant are correct and Friedmann equations are valid. 

The simplest cure (but "embarrassing" from the point of view of the particle physicist), is to add a 
cosmological constant, or dark energy. It should be stressed that this minimal modification of Friedmann 
equations is consistent with all other current cosmological tests and measurements. In the general case, 
the age of the Universe can be related to the expansion history as 

rk) i-to rco 

to= dt= adij= . (59) 

Jo Jo Jo [^ + z)H{z) 

Here we have used Eqs. (l24l) . d27l . d28l . For two components, pressureless matter and dark energy with 
equations of state w, this relation can be written as (see Eq. (l36l ): 

Hoto = / ^ ,^ . (60) 



(1 + z)5/VS^M + ^^De(1 + 



l3ui 



The product Hoto as a function of w, assuming + ^de = 1> is shown in Fig. |9l We see that it 
matches the observational constraints when w — 1 and Vl^.^ w 0.3. 

A discussion of further evidence for dark energy, e.g., related to the problem of the growth of 
density perturbations, can be found in Ref. L32J . 

Redshift - Luminosity Distance relation for Supernovae la. For the two-component energy content 
of the Universe, presureless matter and dark energy, the expression for the luminosity distance, Eq. (I33t . 
takes the form 

^ Ho Jo (l + z')3/V^M + f^DE(l + ^')'" ' 




Fig. 10: Left panel: SN la residual Hubble diagram comparing astrophysical models and models for astrophysical dimming. 
Data (weighted averages in fixed redshift bins) and models are shown relative to an empty Universe model (f2 = 0), adopted 
from Ref. 1111 . Right panel: Joint confidence intervals for (QM,^iA) from SNe la, Ref 1 1 1 1. The dotted contours are the results 
from Ref. 1331 . illustrating the earlier evidence for JIa > 0. The figure is adopted from Ref. 1111 . 

see Eqs. d32b . d33b . d36t . To use this relation as a cosmological test in conjunction with Eq. d34b . one 
has to find a set of standard candles. This is a big challenge in practice, since we have to find very 
bright objects which can be seen from far away. At the same time, all of them should have the same 
luminosity, and we have to be sure that they do not evolve intrinsically. These requirements rule out 
galaxies and quasars. However, supernovae seem to be suitable. They are bright, as bright as the whole 
galaxy at the peak of luminosity, and Type la supernovae appear to be standard candles. These type of 
supernovae are thought to be nuclear explosions of white dwarfs in binary systems. The white dwarf, 
a stellar remnant supported by the degenerate pressure of electrons, accrete matter from a companion 
and its mass increases toward the Chandrasekhar limit of about 1.4 Mq. Near this limit, the degenerate 
electrons become relativistic, which leads to instability and the white dwarf explodes. This physics 
allows the explosions to be calibrated, since instability occur under the same conditions. 

To proceed, I have to remark about the units of flux used by astronomers, which are magnitudes. 
The system is ancient, and has its origin in the logarithmic response of the human eye. The ratio of the 
flux of two objects is then given by a difference in magnitudes; i.e., 

ma - mi = -2.5 logio(i^2M) • (62) 

A smaller magnitude means larger flux. 

Fig. left panel, shows the corresponding SNe la redshift-luminosity distance diagram. Data 
points correspond to magnitudes of SNe la measured at different redshifts. The case of Qm = 1 (red 
curve) is ruled out. The "concordance model" Um = 0.27 and I^a = 0-73 (black dashed curve) is within 
la. For a flat geometry prior, best fit corresponds to i7jv/ = 0.29±q q3 (correspondingly JIa = 0.71). 
Data are inconsistent with a simple model of evolution of SNe la, or dimming due to light absorption by 
dust as an alternatives to dark energy. The shaded area in Fig.^l right panel, corresponds to 68%, 95% 
and 99.7% confidence levels in the (Qm, ^a) plane. 

Constraints from CMB. Tight consti^aints on dark energy, and in a direction in parameter space which 
is "orthogonal" to SNe le constraints, are obtained from fitting the power spectrum of cosmic microwave 
background anisotropics, see the discussion in Section |3 The WMAP data alone rule out the standard 
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Fig. 11: Left panel: Constraints on the equation of state of dark energy in the (w, fiA/) plane for a combination of the CMBR 
+ 2dF + SNe la data sets. Right panel: Constraints on the geometry of the universe in the (^Im, ^a) plane assuming the prior 
h > 0.5. From Ref. E) 



= 1 CDM model by 7a if the prior h > 0.5 is accepted |9|. The resulting confidence levels in 
(^M, ^a) plane are shown in Fig. While the CMBR data alone are compatible with a wide range 
of possible properties for the dark energy, the combination of the WMAP data with either the HST key 
project measurement of Hq, the 2dFGRS measurements of the galaxy power spectrum or the Type la 
supernova measurements requires that the dark energy be Qa = 0.73 it 0.04 of the total density of the 
Universe, and that the equation of state of the dark energy satisfy w < —078 (95% CL) t9J. 



Constraints from gravitational lensing. Gravitational lensing will be discussed in Section |3 Here 
we just note that the analysis of strong lensing of sources with known redshift is sensitive to the value 
of the geometrical cosmological parameters of the Universe. A recent study 135] of the lensing config- 
uration in the cluster Abell 2218 is in agreement with the concordance model. In particular, assuming 
the flat Universe, it gives for the equation of state of dark energy w < —0.85. These constraints are 
consistent with the current constraints derived from CMB anisotropics or supernovae studies, but they 
are completely independent tests, providing nearly orthogonal constraints in the {^m,^a) plane, see 
Fig. EH 



Biggest Blinder - Biggest Surprise. From the point of view of the particle physicis,t the cosmological 
constant just should not be there. Indeed, in quantum field theory, the cosmological constant corresponds 
to vacuum energy, which is infinite and has to be renormalized. 
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cjfc k'^ dk . (63) 



The natural value for the cut-off in this integral is the Plank scale, and then pvac ~ -^^pi ~ 10 GeV^. 
Exact supersymmetry can make this integral vanish. Indeed, in Eq. the contribution of one Bosonic 
degree of freedom is counted. Fermions contribute with an opposite sign, and if there is an equal number 
of Bosons and Fermions with equal masses, the vacuum energy will be zero. However, supersymmetry is 
broken at least at the electroweak scale, and then pvac should not be smaller than ~ ^ 10^ GeV^. 

Before dark energy was detected, it was believed that some yet unknown mechanism reduces the cos- 
mological constant to zero. Zero is a natural number. However, it is hard to understand the smallness 
of the observed value pvac ~ 10"'*^ GeV^. Moreover, there is another pressing issue of fine tuning: 
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Fig. 12: The rotational curves of two sample galaxies. Left panel - M33, adapted from Ref. f361. Right panel - NGC6503, 
adapted from Ref. 1371 . I superimposed with the rotational curves the optical images of corresponding galaxies, approximately 
to scale. 

why the detected value of pvac approximately equals to the energy density of matter at the present epoch 
of cosmological evolution? The ratio of these two contributions scales as and, say, at recombination 
the vacuum energy was only 10^^ of matter energy... Detection of dark energy not only points to a new 
physics, but hints that we are missing SOMETHING very fundamental. 

6. DARK MATTER 

CMBR observations accurately measure the geometry of the Universe, its present expansion rate, its 
composition, and the nature and spectrum of the primordial fluctuations. Nevertheless, the traditional 
cosmological tests are still important. In particular, degeneracies between different parameter sets exist, 
which can produce the same CMBR spectra, and the conclusions drawn do rest upon a number of as- 
sumptions. Below we consider cosmological observations that are independent of the CMB and point to 
the existence of non-baryonic dark matter. 

6.1 DARK MATTER: motivation 

The missing mass is seen on all cosmological scales and reveals itself via: 

• Flat rotational curves in galaxies. 

• Gravitational potential which confines galaxies and hot gas in clusters. 

• Gravitational lenses in clusters. 

• Gravitational potential which allows structure formation from tiny primeval perturbations. 
6.11 Dark Matter in Galaxies 

Galactic rotational curves. Consider a test particle which is orbiting a body of mass M at a distance 
r. Within the frameworks of Newtonian dynamics the velocity of a particle is given by 




(64) 



Outside of the body, the mass does not depend on distance, and the rotational velocity should obey 
the Kepler law, v^ot oc r~^/^. Planets of the Solar system obey this law. However, this is not the case 




Fig. 13: Left panel: the phase space structure of an infall model. Right panel: rotational curves in an infall models. Two curves 
which correspond to different angular momenta are shown. From Ref. 1401 . 



for stars or gas which are orbiting galaxies. Far away from the visible part of a galaxy, rotational curves 
are still rising or remain flat. Two examples ai^e shown in Fig. El An optical image of the M33 galaxy 
is superimposed with its rotational curve, approximately to correct scale. The contribution of visible 
baryons in the form of stars and hot gas can be accounted for, and the expected rotational curve can be 
constructed. The corresponding contributions are shown in Fig. One can see that the data-points 
are far above the contribution of visible matter. The contribution of missing dark mass, which should be 
added to cope with data, is also shown and is indicated as Dark halo. For the rotational velocity to remain 
flat, the mass in the halo should grow with the radius as M(r) oc r, i.e., the density of dark matter in the 
halo should decrease as p(r) oc r~^. 



Halo structure. For direct and indirect dark matter searches it is important to know the phase-space 
structure of the dark halo as well. With dark matter particles that are interactive, a thermal distribu- 
tion over velocities would eventually be established. However, in conventional cold dark matter models, 
particles are non-interacting, except gravitationally. Binary gravitational interactions are negligible for 
elementary particles, and the resulting phase-space distributions are not unique, even for stationary equi- 
librium states, and even if flat rotational curves are reproduced. 

1. The simplest self-gravitating stationary solution which gives flat rotational curves corresponds 
to an "isothermal sphere" with Maxwellian distribution of particles over velocities: 

n(f,r;) = n(r) e-'^'/^o . (65) 

Solution of the equation of hydrostatic equilibrium can be approximated by the density profile 

P(^) = /-I ^° 9x > where x = r/rc . (66) 
(1 + x^j 

It should be stressed that the distribution Eq. (I65t . in contrast to a distribution in real thermal equilibrium, 
depends on particle velocities, not on their energies. Such distributions may arise in time-dependent 
gravitational potential as a result of coUisionless relaxation. 

2. There exist several density profiles which are empirical fits to numerical simulations, e.g. 
Navarro, Frenk & White (NFW) profile OH and Moore et al. profile EH 

P(0 = ^ .2 NFW, (67) 



X (1 + x)^ 
Po 

3.3/2 (1 + ^3/2) 



p{r) = ^° — — Moore et al. (68) 




Fig. 14: Left panel: one of the detected microlensing effects, see Ref. 1441 . Right panel: schematic view of gravitational lensing 
by point mass. 



4. In the cold dark matter model, the distribution of particles in the phase space during initial 
linear stage prior to structure formation corresponds to thin hypersurface, v = Hr (or line in the Hubble 
diagram). Since during coUisionless evolution the phase-space density conserves, at the non-linear stage 
the distribution will still be a thin hypersurface. If can be deformed in a complicated way and wrapped 
around, but it cannot tear apart, intersect its own folds, puff up or dissolve. The corresponding phase- 
space distribution for the case of spherical symmetry is shown in Fig.^J left panel. With time, non-linear 
structure grows, and the infall of new particles continues. This manifests itself as a growth of turnaround 
radius (which is a surface where v = 0; the turnaround radius of our Galaxy is at 1 Mpc, see Ref. |41 1) 
and as an increasing number of folds inside turnaround. The energy spectrum of dark matter particles 
at a fixed position will be discrete, see Fig. left panel, where several velocity peaks are indicated at 
intersections of the vertical dashed line, r = const, with phase-space sheets. The overall shape of the 
spectrum also changes compared to an isothermal distribution. This may be important for direct dark 
matter searches. 

The infall model reproduces flat rotational curves, see Fig.^J right panel. There is one interesting 
difference, though; rotational curves of the infall model have several small ripples in the region where 
the curve is flat. These ripples appear near the surfaces where ?; = 0. In principle, they may be detectable 
||42 43 1 and then it will give a clear, unique signature of the presence of dark matter in galactic halo, (as 
opposed to models which try to explain apparent violation of Kepler's law by modification of gravity). 

Existence of such a folded structure is a topological statement. However, in the inner halo the 
number of folds is very large, and limited resolution makes the distribution indistinguishable from, say, 
isothermal. It is not clear at which distances the description of halo in terms of the infall becomes 
appropriate. But for sufficiently isolated galaxy, in regions closer to the outer rim of the halo, where the 
number of folds is still small, signatures of the infall should exist, and they do exist in our Galaxy 1411 . 

Baryonic Halo Dark Matter ? No. Already CMBR alone tells us that there should be non-baryonic 
dark matter, see Table I. BEN and CMB agree on = 0.04, however, contribution of stars amounts 
only to Ostars = 0.005. There should be dark baryons hiding somewhere. Can it be that the whole, or 
at least some part, of the halo dark matter is comprised of dark baryons in the form of non-luminous 
objects? Candidates are Jupiter like planets, brown dwarfs (which are undersized stars, too light to ignite 
thermonuclear reactions), or already dead stars (white dwarfs, neutron stars, or even black holes). This 
class of objects got the acronym MACHO, from MAssive Compact Halo Objects. Special techniques 
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Fig. 15: 95% confidence limits on the MACHO fraction of the standard local halo density. Green and blue lines - results of 
EROS |45 1 and MACHO [46| microlensing collaborations. Red line - constraint from the absence of distortion of distribution 
of binary stars in angular separation 147 J . Note that the microlensing exclusion curve extends outside of the plotted range and 
up to M ~ 1O~^M0, see Ref. jlU. 



based on gravitational lensing were developed for MACHO searches. These searches were successful, 
but by now it is clear that MACHOs cannot comprise the whole dark matter, as their fraction of DM halo 
is restricted to be < 50%. Since MACHOS are the only type of dark matter which has been detected, let 
us consider the issue in some more detail. 

Consider the light deflection by a point mass M. If the impact parameter ^ is much larger than the 
Schwarzschild radius of the lens, ^ » 2GM, then General Relativity predicts that the deflection angle 
of a light ray, a, is 

AGM 

a = —— . (69) 

This is twice the value obtained in Newtonian gravity. If the lens happens to be on the line which connects 
the observer and a source, the image appears as a ring with the radius (Einstein ring radius) 

4 = 4GM^-^, (70) 

see Fig. [2] right panel. If the deflector is displaced from the line of sight by the distance r, then instead 
of the ring, an odd number of images will appear. If the images cannot be observed separately, because 
the resolution power of the telescope is not sufficient, then the only effect will be an apparent brightening 
of the source, an effect known as gravitational microlensing. The amplification factor is 

■ (7.) 



where u = r/vE- If the lens is moving, the distance r will be changing with time, and the image of 
the background star will brighten during the closest approach to the line of sight. If the galactic halo is 
filled with MACHOs, this may happen occasionally for some of background stars. The typical duration 
of the light curve is the time it takes a MACHO to cross an Einstein radius. At ~ ve/v, where v ~ 10^'^ 
is typical velocity in the halo. If the deflector mass is IMq, the average microlensing time will be 3 
months, for 10"^ Mq it is 9 days, for 10~^ Mq it is 1 day, and for 10"*^ M© it is 2 hours. 

An optical depth for microlensing of the galactic halo is approximately r 10~^. Thus, if one 
looks simultaneously at several millions of stars during an extended period of time, one has a good 
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Fig. 16: X-ray (left panel) and optical (middle panel) images of the Abell 2029 cluster of galaxies. Right panel: ratio of total 
enclosed cluster mass to light in A2029, from Ref. |49l. 



chance of seeing at least a few of them brightened by a dark halo object. The first microlensing events 
were reported in 1993. Nowadays, there are more than half a hundred registered events. One of them is 
shown in Fig. [H] However, derived optical depth is not sufficient to account for all dark matter in the 
Galaxy halo. 95% confidence limits on the MACHO fraction of the standard local halo density is shown 
in Fig. [B] 

Since MACHOs cannot account for the mass of the dark halo, non-baryonic dark matter should be 
present out there. 



6.12 Dark Matter in Clusters of Galaxies 

Already in 1933, F. Zwicky BHll deduced the existence of dark matter in the Coma cluster of galaxies. 
Nowadays, there are several ways to estimate masses of clusters, based on the kinetic motion of member 
galaxies, on X-ray data, and on gravitational lensing. These methods are different and independent. In 
the dynamical method, it is assumed that clusters are in virial equilibrium, and the virialized mass is 
easily computed from the velocity dispersion. In X-ray imaging of hot intracluster gas, mass estimates 
are obtained assuming hydrostatic equilibrium. Mass estimates based on lensing are free of any such 
assumptions. All methods give results which are consistent with each other, and tell that the mass of the 
luminous matter in clusters is much smaller than the total mass. 



Kinetic mass estimates. Those are based on the virial theorem, (-Epot) + 2(£'kin) = 0. Here (-Ekin) = 
N{mv'^)/2 is averaged kinetic energy of a gravitationally bound object (e.g. cluster of N galaxies) and 
(Epot) = —N'^{m?) /2{r) is its averaged potential energy. Measuring the velocity dispersion of galaxies 
in the clusters and its geometrical size gives an estimate of the total mass, M = N{m). 

M ~ ^ . ,72) 

The result can be expressed as mass-to-light ratio, M/L, using the Solar value of this parameter. For the 
Coma cluster, which consists of about 1000 galaxies, Zwicky 1.481 has found 

^^ 300/.^. (73) 

Modern techniques end up with very much the same answer. M/L ratios measured in Solar units in 
central regions of galaxies range from a few to 10 in spirals and large ellipticals. If clusters are large 
enough systems for their M/L to be representative of the entire Universe, one finds i50l 



Qm ~ 0.2 - 0.3 



(74) 
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Fig. 17: Left panel: an image of the cluster Abell 2218 taken with the Hubble space telescope (see Ref. [52]). Spectacular arcs 
resulting from strong lensing of background galaxies are clearly seen. 

Mass estimates based on X-rays. Mass is also traced in clusters of galaxies by hot gas which is visible 
in X-rays. Assume hot gas is in thermal equilibrium in a gravitational well created by a cluster. Then its 
density distribution Pg{r) and pressure Pg{r) satisfy 



Observationally, the gas density follows from the X-ray luminosity. Gas temperature can be measured 
from the shape of X-ray spectrum. By measuring the temperature profile of a gas, one can reconstruct 
the gas pressure Pg{r). In this way, the radial run of mass can be deduced. 

For example, detailed modeling l49l of Abell 2029, which is shown in Fig.^l leads to the conclu- 
sion that the cluster is dark matter dominated all the way into its core. After subtracting the contributions 
of stars and hot gas into the mass budget, the density profile of dark matter can be reconstructed. It 
agrees with NFW dark matter profile, Eq. d68l . poc \/x{l + x^), where x = r/vg and = 540 kpc. 
The agreement is remarkably good on all scales measured, 3 - 260 kpc. Baryons contribute fh 14% 
into the total mass of the cluster. Assuming universal baryon mass fraction and Q;, from big bang nucle- 
osynthesis, this also gives $7^ ~ ^-m/ fb ~ 0.29 for the total mass budget in the Universe, in agreement 
with other current estimates. 

The same methods can be employed for studies of dark matter in large elliptical galaxies. In 
Ref. II5TI the mass profile of the elliptical galaxy NGC 4636, based on the temperature of hot interstellar 
gas, was obtained for distances from 0.7 to 35 kpc. It was found that the total mass increases as radius to 
the power 1.2 over this range in radii, attaining a mass-to-light ratio of 40 solar masses per solar visual 
luminosity at 35 kpc. As much as 80% of the mass within the optical half-light radius is non-luminous 
in this galaxy. 

Gravitational Lensing. Gravitational lensing allows direct mass measurement without any assump- 
tions about the dynamical state of the cluster. The method relies on the measurement of the distortions 
that lensing induces in the images of background galaxies. As photons travel from a background galaxy 
to the observer, their trajectories are bent by mass distributions, see Fig. \n\ right panel. Consider the 
deflection by a point mass M. For impact parameter ^ which is much larger than the Schwarzschild 
radius of the lens, ^ 2GM, the deflection angle a is given by Eq. (l69t . If the gravitational field is 
weak, the deflection angle of an ensemble of point masses will be the vectorial sum of the deflections 
due to individual lenses. 

A reconstruction of lens geometry provides a map of the mass distribution in the deflector. For a 
review of the method see e.g. Ref l53l . The images of extended sources are deformed by the gravitational 
field. In some cases, the distortion is strong enough to be recognized as arcs produced by galaxy clusters 
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Fig. 18: The constraints on cosmological parameters obtained from the study of the lensing configuration in the Abell 2218 
cluster of galaxies 1 35 1. 



serving as a lens, see Fig.EJ left panel. For the cluster A 2218, shown in this figure. Squires et al. l54l 
compared the mass profiles derived from weak lensing data and the X-ray emission. The reconstructed 
mass map qualitatively agrees with the optical and X-ray light distributions. A mass-to-light ratio of 
M/L = (440 ± 80) /i in solar units was found. Within the error bars the radial mass profile agrees with 
the mass distribution obtained from the X-ray analysis, with a slight indication that at large radii the 
lensing mass is larger than the mass inferred from X-rays. The gas to total mass ratio was found to be 
Mgas/Mtot = (0.04 ± 0.02) /1-3/2. 

Interestingly, the analysis of multiple images of several sources with known (and significantly 
different) redshift produced by a cluster lens is sensitive to the value of the geometrical cosmological 
parameters of the Universe. Study |35 | of the lensing configuration in the cluster Abell 2218 gives 
< Q.M < 0.30 assuming a flat Universe, and < Q.m < 0.33 and w < —0.85 assuming fiat 
Universe with dark energy, see Fig. \W\ These constraints are consistent with the current constraints 
derived with CMB anisotropics or supernovae studies, however this method is a completely independent 
test, providing nearly orthogonal constraints in the (Om, ^^a) plane. 

6.2 Structure formation and DM 

By present time the structure in the Universe (galaxies and clusters) is formed already, the perturbations 
in matter 6p/p > 1. However, the initial perturbations were small Sp/p 10^^. Perturbations 
do not grow in the radiation dominated epoch, they can start growing only during matter domination 
6p/p ~ a = 1/z. Moreover, baryonic plasma is tightly coupled to radiation, therefore perturbations 
in baryonic matter start to grow only after recombination. For the same reason, initial perturbations in 
baryons at the time of recombination equal to fluctuations in CMBR. If baryons were to constitute the 
only matter content, then perturbations in matter at present time would be equal to 

today = ^rec rec ~ 10 , (76) 

p p 

where z^ec ~ 1100 is the redshift of recombination. This is one of the strongest and simplest arguments in 
favour of non-baryonic dark matter. Structure has had time to develop only because perturbations in non- 
baryonic dark matter have started their growth prior to recombination. Baryonic matter then "catch up" 
simply by falling into already existing gravitational wells. If one aims to explain things by modification 
of gravity, one has to explain not only the flat rotational curves in galaxies and the presence of dark 
matter in galaxy clusters, but has also to provide the accelerated growth of structure from recombination 
till present, in a consistent way. 
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Fig. 19: Left panel: The constraints on axion parameters. Red blocks - various cosmological and astrophysical constraints; 
yellow blocks - exclusion regions obtained in dedicated dark matter search experiments; green blocks - the allowed regions 
in two cosmological scenarios. Right panel: Solid curve - the exclusion limits on the coherent WIMP-nucleon scalar cross- 
section obtained by CDMS collaboration in the direct dark matter search experiment; the parameter space above the curve is 
excluded at the 90% C.L. These limits are inconsistent with the DAMA 3a signal region |73| (closed dashed contour) if it is 
due to scalar coherent WIMP interactions. Also shown are limits from CDMS at SUF (dotted line). The typical predictions of 
supersymmetric models are shown in yellow. Adapted from Ref. 1741 . 



6.3 Non-baryonic Dark Matter Candidates 

There is no lack for dark matter candidates in particle physics models. Some of them appear naturally, 
and were motivated by some other reasoning, not related to the dark matter problem. They are the leading 
candidates and are listed below. 

• Axion. Has a mass m ~ 10~^ eV. Appeared 1551 B6l as a by-product of a suggested solution 
of the strong CP problem via a global C/(l) Peccei-Quinn symmetry |57|. The axion picks up a 
small mass in a way similar to the pion when chiral symmetry is broken. The parameters of these 
two particles are related, in particular rua ~ rn.„-f.n-/ fa, where fa, is the axion decay constant, and 
determines also the strength of the axion coupling to other particles. There are tight astrophysical 
and cosmological bounds on fa which leave only a narrow window, 10^*^ GeV < /a ^ 10^^ GeV, 
see Fig.^1 the left panel. For the review of axion physics and searches see e.g. Refs. I58ll59ll . 

• Neutrino, m ~ 0.1 eV. The only dark matter candidate which is known to exist. For this reason 
we discuss the neutrino in some more detail below. For the review of neutrino cosmology see e.g. 
Ref. [60 1, and for the neutrino astrophysics see e.g. Ref. 16X1 . While the neutrino is cosmologically 
important, it cannot resolve the dark matter problem. 

• Mirror matter, M ~ 1 GeV. Does not belong to the list of the most popular candidates, but is 
attractive as an example of a model l62l l63l l64l where the approximate equality of baryonic and 
non-baryonic contributions into the energy balance of the Universe is attempted to be explained 
naturally, and not as a result of fine-tuning of model parameters. 

• WIMP, m ~ 100 GeV. The most popular candidate, a natural outcome of supersymmetry. The 
lightest supersymmetric particle (or LSP) is naturally stable and would have interesting cosmo- 
logical abundance. Known also under the names Neutrallino (dark matter has to be color and 
electrically neutral) and WIMP (from Weakly Interacting Massive Particle). For recent reviews 
see e.g. Refs. li65n66J . The current status of direct and indirect WIMP searches is reviewed in 



Ref. ll67ll . The new limits obtained by the CDMS collaboration and not reflected in cited reviews 

are shown in Fig. ^1 the right panel. 
• WIMPZILLA, m ~ 10^'^ GeV. The newcomer, was initially motivated as a solution of the Greizen- 

Zatsepin-Kuzmin puzzle of Ultra-high energy cosmic rays f 68ll69l . The popularity was boosted 

by the observation that cosmologically interesting abundance is created naturally, just as a sole 

consequence of Universe expansion [70 71 1. 
Non-baryonic dark matter model building and searches makes and extensive subject on its own. There 
are many dedicated excellent reviews, I cannot list them all, see e.g. some earlier |72| and the latest 
I65ll66l l67 1 one. For this reason, and because of space limitations, I will not describe non-baryonic dark 
matter in all its variety, instead I'll spend some time on simple and universal relations. 



Cosmological density of neutrino. Here we calculate the abundance of particles which were once in 
thermal equilibrium with the rest of cosmological plasma. Let us first consider the case of the neutrino. 

Comparing the weak interaction rate, F ~ T^/M^, to the expansion rate, H ~ T^/Mpi, one 
finds that neutrino are in thermal equilibrium at temperatures T » 1 MeV and decouple from the rest of 
plasma at lower temperatures. (One can do this in full detail, see e.g. l60l .') Therefore, standard model 
neutrinos, which have small masses, decouple when they are still relativistic. The number density of 
neutrino at this time is given by Eq. (fT4b . Below this temperature, neutrinos are no longer in thermal 
equilibrium with the rest of the plasma, and their temperature simply decreases as T oc 1/a. However, 
the cosmological background of photons is heated up by the e+e^ annihilations. Let us find a relation 
between T^, and Tj, which will also give the relation between n,^ and n^. 

The annihilation reaction rate is much faster than the expansion of the universe, therefore this 
process is adiabatic, and entropy in comoving volume is conserved, = const, see Eq. (fm) . Before 

annihilation g* = + ge ■ (7/8) = 2 + 4 • (7/8) = 11/2. After annihilation ff* = fii-y = 2. Since before 
annihilation T^y = T^, we find T^, = (4/11)^^^ T^, and for one neutrino flavour we have 

3 

riu = — n^ ■ (77) 

Here we used Eq. (fT4l) and g^ = gu = 2, since right handed neutrino do not contribute (are not excited) 
even if neutrino has a small Dirac mass, see Ref. |60|. As a consequence of e+e" annihilation neutrino 
temperature is lower. With = 2.728 K we find = 1.947 K and = 115 cm'^. 

At temperatures larger than neutrino mass, Tjy » m^, in the standard model, assuming no chemi- 
cal potential, we find 

^u = 3 G) (n) ^7^0-68 5^7 ' (78) 

where the factor of 3 corresponds to the number of neutrino flavours. This result allows us to find the 
epoch of equal matter and radiation densities 

Assume that by now neutrino became non-relativistic, i.e. their masses are larger than the present 
temperature. In this case, neutrino energy density is given by p^, = niui riyi. Since it has to be smaller 
than Q.rn Pc, we have the constraint f75l 

rriyi < 94 VL^h^ eV = 12 eV . (80) 

i 

For dark matter particles to boost the structure formation, their typical velocities squared at the time 
of recombination should be smaller than the depth of typical gravitational wells, <^ 10^^. In other 
words, the dark matter should be cold. This is not the case for particles as light as those which satisfy 
the bound Eq. dSOl . Neutrino can make up dark matter, but it will be hot dark matter. 



Neutrino mass is pinned down. Free streaming of relativistic neutrinos suppresses the growth of fluc- 
tuations until V becomes nonrelativistic at z ~ mj/3To ^ 1000 {rrij/eY). This effect of free- 
streaming is not seen in the data and therefore only small corrections due to light neutrino are allowed in 
the standard CDM picture. Combined CMBR and LSS analysis yields the constraint ||9j 

^uh^ = < 0.0076 , (81) 

93.5 eV 

which translates into the upper bound 

< 0.7 eV (95% CL) . (82) 

i 

On the other hand, atmospheric neutrino oscillations provide a lower bound on the heaviest neutrino 
mass, since \/6rn^^ ~ 0.03 eV. Combining these two limits 

0.03 eV < mheaviest < 0.24 eV (83) 
we see that the heaviest neutrino mass is now known to within an order of magnitude fJ6\. 



Can neutrino make up a galaxy iialo? By z 1 the neutrino quanta satisfying the mass bound 
Eq. (l83b became sufficiently non-relativistic to make their way into gravitational wells. The question 
arises, can neutrino at least make up the dark halos and be responsible for flat rotational curves? The 
answer to this question is: no. To prove it, let us assume that neutrino does build up a dark matter halo 
with a flat rotational curve 

^DM - • C«4) 

We can express the energy density pjy^^ = through the integral of phase space density over the 
momenta 



p, = j^ J d-'kn{k,r) , (85) 

But for fermions, the phase-space density, n{k, r) = n{E), should obey the Pauli exclusion principle, 
n{E) < 1. Combining Eqs. (l84l i and (l85l . we find m'^v'^ ~ Mp^v^ /r"^, or 

/lOO kms-i\^/^ /lkpc\^/^ 

rriu > 120 eV — — (86) 

V Vrot J \ rc J 



For dwarf galaxies this constraint (the Tremaine-Gunn limit iTTTI ) reads > 500 eV, and we arrive to 
contradiction with Eq. dSOl . which becomes even stronger when compared with Eq. (l83l . 



Cosmological density of otiier thermal relics. Assume now some weakly interacting particle has a 
mass larger 1 MeV and decouples when it is non-relativistic. The equilibrium number density will be 
Boltzmann-suppressed in this case, by the exponent exp(— m/T). The weak interaction cross-section 
implies a ~ m^/M^, if m <C Mw- Repeating calculations for abundances, one finds that in this case 
Qmh'^ 3 (1 GeV/m)^, i.e. a correct cosmological abundance of dark matter would be achieved for 
m ~ 5 GeV. 

On the other hand, if m » Mw, the annihilation crosssection becomes a ~ 1/m^ and one finds 
Qmh'^ ~ (m/1 TeV)^, i.e. the correct cosmological abundance of dark matter is achieved for m w 300 
GeV. Using field-theoretical unitarity and the observed density of the Universe, it can be shown that a 
stable elementary particle which was once in thermal equilibrium cannot have a mass greater than 340 
TeV |,78J. 



Cosmological density of non-thermal relics. The mass of non-thermal rehcs can be much larger than 
0(10^ ) TeV without violating unitarity bound; it can also be much smaller than than 0(1) GeV and dark 
matter still will be cold, as required by observations. 

1. Axions. Very light scalar particles, like axions, are created in a state of coherent oscillations. 
This can be viewed also as a Bose-condensate. To illustrate the general idea, let us consider a scalar field 
with potential V{(j)) = rn^cj)^ /2. The equations of motion for the Fourier modes with a momentum k in 
an expanding Universe are 

'<Pk + 3H^k + {k^ + m^)(kk = . (87) 

Since the term oc H can be understood as friction, amplitude of those modes for which QH^ ^> {k'^+m?') 
(almost) does not change with time. The oscillations of modes with a given k commence when H 
becomes sufficiently small, QH^ <^ (fc^ + m^). Oscillating modes behave like particles, and their 
amplitude decreases with expansion. Since modes with the largest k start oscillations first, they will have 
the smallest amplitude and the field becomes homogeneous on a current horizon scale. Modes with all 
k will be already oscillating when 3H m, and will behave like cold dark matter since then. Note 
that the field will be homogeneous on the horizon scale at this time, but may be inhomogeneous on 
larger scales. This may lead to formation of dense clumps, "axion mini-clusters" [19. 80 1 of the mass 
M ~ 10~^^ Mq IsTTI . Note also that in the case of axions, one has to take into account the dependence 
of m on temperature T. Solving 3H{T) = m{T) one finds J7axion ~ 1 when fa ~ lO^^QgV I82U83II . 

2. Superheavy dark matter. Non-conformal quantum fields cannot be kept in a vacuum in an 
expanding universe. This can be understood on the example of a scalar field, Eq. dSTl . In conformal 
time, Eq. (l24l) . and for rescaled field, Uk = (pk cl, the mode equations take form of an oscillator equation 

ilk + oJkUk = , (88) 

with time-dependent frequency 

u;l = k'^ + a'^m'^ --(1-60 ■ (89) 
a 

The constant ^ describes the coupling to the scalar curvature, corresponding term in the Lagrangian is 
^R(j). The case of ^ = corresponds to minimal coupling (Eq. dSTl was written for this case), while 
^ = 1/6 is the case of conformal coupling. Equations for massless, conformally coupled quanta are 
reduced to the equation of motion in Minkowski space-time. Particle creation does not occur in this 
case. For massive particles, conformal invariance is broken and particles are created regardless of the 
value of ^. Let us consider the case of ^ = 1/6 (the general situation is considered in |84|). It is the 
particle mass which couples the system to the background expansion and serves as the source of particle 
creation in this case. Therefore, we expect that the number of created particles in comoving volume is 
oc and the effect is strongest for the heaviest particles. In fact, stable particles with m > 10^ GeV 
would overdose the Universe in the standard "pre-inflationary" Friedmann model (tT\. Inflation cuts the 
particle production and QsR ~ 1 if ?n. > 10^^ GeV and reheating temperature is T ~ 10^ GeV, which is 
the value of reheating temperature compatible with supergravity models Il70ll7ni84l . 

7. BASICS OF INFLATION 

In frameworks of "classical" cosmology and assuming no fine-tuning, one concludes that a typical uni- 
verse should have had Plankian size, live Plankian time and contain just a few particles. This con- 
clusion is based on the observation that Fridmann equations contain a single dimension-full parameter 
Mpi ~ 10^^ GeV, while dimensionless parameters naturally ai^e expected to be of order unity. Yet, the 
observable Universe contains 10^'^ particles in it and had survived 10^^ Plankian times. Where does it 
all came from? In other words, why is the Universe so big, flat (Qq « 1) and old (t > 10^^ years), 
homogeneous and isotropic {6T/T ~ 10^^), why does it contain so much entropy (S > 10^^) and does 



not contain unwanted relics like magnetic monopoles? These puzzles of classical cosmology were solved 
with the invention of Inflation [85 86||W1|88||^. All these questions are related to the initial conditions 
and one can simply postulate them. The beauty of Inflation is that it prepares these unnatural initial con- 
ditions of Big Bang, while the pre-existing state (which can be arbitrary to a large extent) is forgotten. 
Inflationary theory came with unplanned bonuses. Not only does the Universe become clean and homo- 
geneous during inflation, but also the tiny perturbations necessary for the genesis of galaxies are created 
with the correct magnitude and spectrum. Below we consider the basics of inflationary cosmology. 

7.1 Big Bang puzzles and Inflationary solutions 

Horizon problem and the solution. The size of a causally connected region (horizon) scales in pro- 
portion to time, i?H oc t. On the other hand, the physical size of a given patch grows in proportion to the 
scale factor, Rp oc a{t) (xf^. The exponent 7 depends upon the equation of state, 7 = 1/2 for radiation 
and 7 = 2/3 for matter dominated expansion. In any case, for the "classical" Friedmann Universe 7 < 1 
and the horizon expands faster than volume. Take the largest visible patch today. It follows that in the 
past its physical size should have been larger than the horizon size at the time (since they are equal today) 
and therefore this patch should have contained many casually disconnected regions. For example, as we 
have found in Section|3 the angular size of horizon at the moment of last scattering is « 2°, see Eq. (I55t . 
which tells us that we observe 10^ causally disconnected regions at the surface of last scattering. The 
question arises, why is the Universe so homogeneous at large scales? 

This problem can be solved if during some period of time the volume had expanded faster than the 
horizon. During such a period, the whole visible Universe can be inflated from one ("small") causally 
connected region. Clearly, this happens if 7 > 1, which means a > 0. Either of these two conditions 
can be used as definition of an inflationary regime. Using the Friedmann equation Q we find that the 
inflationary stage is realized when p < —p/3. In particular, if p = —p the energy density remains 
constant during expansion in accord with the first law of thermodynamics, Eq. (Q, and the physical 
volume expands exponentially after a few Hubble expansion times, a{t) = e^*, see Eq. 

Curvature problem and the solution. The Friedmann equation ^ can be re-written as 



Here we immediately see the problem: during matter or radiation dominated stages, decreases (in 
general, this happens for any expansion stage with a < 0), therefore J7 is driven away from unity. To 
observe Q 1 today, the observer has to live in a universe with extreme initial fine-tuning, say at the 
epoch of nucleosynthesis, when temperature was T ~ 1 MeV, one finds \n{t^s) - 1| < 10^^^ and even 
stronger tuning is required at earlier epochs. A possible solution is obvious: accelerated expansion d > 
increases d and therefore drives Q{t) to unity. A robust, crucial and testable prediction of inflationary 
cosmology is a flat Universe, Q, = 1. 

The problem of Entropy. As we know already, the energy of a vacuum, p = —p, stays constant 
despite the expansion. In this way, room for matter full of energy could have been created. The vacuum 
energy is converted into particles and radiation at some later stage and, in particular, the observed huge 
entropy is created. Potentially, this mechanism works for any inflationary scenario, since the product pa"^ 
is guaranteed to grow whenever d > 0. However, the important question is whether a graceful exit out 
of the inflationary stage and successful reheating is possible. In practice, this issue has killed a number 
of inflationary models. Remarkably, the original model by A. Guth L86J had being ruled out precisely on 
these grounds 8901. 

Inflation has to continue for a sufficiently long time for the problems of horizon, curvature and 
entropy to be solved. All these give roughly the same condition on the number of required "e-foldings" 




(90) 



of inflation l86l and we consider here a (simplified) derivation based on entropy. A precise condition can 
be found e.g. in Ref. |9T]| . Multiplying the current temperature in the universe by its visible size we find 
Taxo ~ 10^°, where xo is the comoving size of the present horizon. The product Ta conserves (up to 
the change in the number of relativistic degrees of freedom, which we neglect here) since the Universe 
expansion is adiabatic after the end of Inflation, see Eq. (IIH . Let Tr denote the reheating temperature 
and = aj/oj the number of inflationary "e-foldings", where is the value of scale factor at the end 
of inflation and Oj at its beginning, respectively. We also want at least the whole visible universe to be 
inflated out of a single causally connected patch, which gives aiXo ~ where Hi is the value of the 
Hubble parameter during Inflation. All this gives the condition^ 

^ > 10^0 . (91) 

In popular models of Inflation the ratio Tr/Hi is within a couple orders of magnitude from unity, and we 
find N Z 65. 



7.2 Models of Inflation 

Consider T^j^ for a scalar field 93 

Tfj,u = d^if duif - g^u (92) 

with the Lagrangian : 

L = d^^d^^-V{^) . (93) 

In a state when all derivatives of are zero, the stress-energy tensor of a scalar field simplifies to = 
V{(p) 5^. This corresponds to a vacuum state. Indeed, comparing with Eq. Q, we find V = p = —p. 
There are two basic ways to arrange (p const and hence to imitate the vacuum-like state. 

1. Consider the potential V{(j)), which has a local minimum with a non-zero energy density sep- 
arated from the true ground state by a potential barrier |86|. A universe which happened to be trapped 
in the meta-stable minimum will stay there for a while (since such a state can decay only via subbar- 
rier tunneling) and expansion of the universe will diminish all field gradients. Then the Universe enters 
a vacuum state. This model is ruled out since the inhomogeneities created during the phase transition 
which terminates the inflationary phase are too large [90|. However, the model is good for illustration 
purposes. The frequently asked question is: how can it be that the energy density stays constant despite 
the expansion? In the model with local potential minimum the energy cannot decrease (classically) below 
the local minimum value, and therefore it has to stay constant despite the expansion. 

2. A. Linde was first to realize that things work in the simplest possible setup fE^\. Consider the 
potential 

V{4>) = ^m2</,2 . (94) 

The equation of field motion in an expanding Universe is (/)+3ff(/)+m^(/) = .If H ^ m, the "friction" 
is too big and the field (almost) does not move. Therefore time derivatives in T^^, can be neglected, and 
inflation starts (in a sufficiently homogeneous patch of the Universe). A Hubble parameter in this case 
is determined by the potential energy, H m(j)/Mpi, and we see that inflation starts if the initial field 
value happens to satisfy (/) > Afpi- During inflationary stage the field slowly rolls down the potential 
hill. This motion is very important in the theory of structure creation; inflation ends when cf) Mpi. At 
this time, field oscillations start around the potential minimum and later decay into radiation. In this way 
matter was likely created in our Universe. 

'strictly speaking in this relation Tr is not the real temperature in a state of thermal equilibrium, but T,. ~ pl^*, where p,. 
is the energy density at the moment when the expansion becomes dominated by relativistic particles. 



7.3 Unified theory of Creation 

During Inflation and by its end the Universe was in a vacuum-like state. We have to figure out how this 
"vacuum" was turned into the matter we observe around us, and how primordial fluctuations which gave 
rise to galaxies were created. Fortunately, these problems can be formalized in a nice and unified way. 
Basically, everything reduces to a problem of particle creation in a time-dependent classical background. 
On top of every "vacuum" there are fluctuations of all quantum fields which are present in a given 
model. This bath of virtual quanta is indestructible, and even Inflation cannot get rid of it. Being small, 
fluctuations of any field obey an oscillator equation 

Uk + [k'^ + ml^iri)] Ufc = 0, (95) 

here are amplitudes of fluctuating fields in Fourier space. Effective mass becomes time dependent 
through the coupling to time-dependent background. Because mcff is time dependent, it is not possible 
to keep fluctuations in a vacuum. If one arranges to put oscillators with momentum k into the vacuum at 
one time, they will not be in vacuum at a latter time because positive and negative frequency solutions 
mix, see below. Several remarks are in order. 

• Eq. (l95l is valid for all particle species. 

• The equation looks that simple in a conformal reference frame ds^ = a{rj)'^ {d-rf — dx^). (Every- 
where in this chapter a "dot" means derivative with respect to 77.) 

• Of particular interest are ripples of space-time itself: curvature fluctuations (scalar fluctuations of 
the metric) and gravity waves (tensor fluctuations of the metric). 

• Effective mass mcff can be non-zero even for massless fields. Gravitational waves give the simplest 
example |92|, with m^^ = —a/ a. The effective mass for curvature fluctuations has a similar 
structure m^g = —z/z, but with a being replaced by 2; = a(j)/H, see Refs. Il93ll94ll95ll96l . 

• For conformally coupled, but massive scalar rricg = mo a{v)- 

Note that creation was only possible because nature is not conformally-invariant. Otherwise, rrics is time- 
independent and vacuum remains vacuum forever. There are two important instances of time varying 
classical background in cosmology: 

• Expansion of space-time, a{7]). 

• Motion of the inflaton field, (/){r]). 

Both can be operational at any epoch of creation: 

• During inflation (superhorizon size perturbations). 

• While the inflaton oscillates (reheating). 

During inflation superhorizon size perturbations of metric are created, which give seeds for Large 
Scale Structure (LSS) formation and eventually lead to the formation of galaxies, and therefore of the 
Solar system and all the rest which we can see around us. During reheating matter itself is created. 
Overall, there are four different situations (two sources times two epochs). If coupling to the inflaton is 
not essential, the corresponding process will be called "pure gravitational creation" in what follows. 

There are several primary observables which can be calculated out of and further used for 
calculation of quantities of interest. Most useful are: 

• The particle occupation numbers, n^. Integration over d'^k gives the particle number density. 

• The power spectrum of field fluctuations, = u^^Uk- Integration over d'^k gives the field variance. 
Depending on physical situation, only one or the other may have sense. The particle number in a comov- 
ing volume is useful because it 

• is adiabatic invariant on sub-horizon scales (or when m > H); 

• allows to calculate abundances of various relics, e.g. dark matter. 

But it has no meaning at super-horizon scales when m < H. The power spectrum and/or field variance 
is useful because it 



• does not evolve on super-horizon scales if m < H; 

• allows to calculate density perturbations generated during inflation; 

• is crucial for dynamics of phase transitions; 

• helps to calculate back-reaction in a simple way (Hartree approximation). 
But Pfc evolves on sub-horizon scales and when m > H. 

Let me start with the discussion of metric perturbations. 



Gravitational creation of metric perturbations. As an important and simple example, let us consider 
quantum fluctuations of a real scalar field, which we denote as if. It is appropriate to rescale the field 
values by the scale factor, ip = (p/a{7]). This brings the equations of motion for the field (j) into a simple 
form of Eq. ( l95l . As usual, we decompose (j) over creation and annihilation operators 6^ and 



(96) 

Mode functions satisfy Eq. (|95J). In what follows we will assume that ip is the inflaton field of the 
"chaotic" inflationary model, Eq. ( |94T ). During inflation H ^ m and H const. So, to start with, we 
can assume that is a massless field on the constant deSitter background. (The massive case can be 
treated similarly, but analytical expressions are somewhat more complicated and do not change the result 
in a significant way. Corrections due to change of H can also be taken into account, and we do that later 
for the purpose of comparison with observations.) With a constant Hubble parameter during inflation the 
solution of Friedmann equations in conformal time is 

a{v) = (97) 

and the equation for mode functions of a massless, conformally coupled to gravity = 0), scalar field 
takes the form 

2 2 

iik + k Uk ^ Ufc = . (98) 

Solutions which start as vacuum fluctuations in the past {rj —>■ — oo) are given by 

Uk = ^= 1 ± ^ • (99) 
V2k V kr]J 

Indeed, at r] — oo the second term in the parentheses can be neglected and we have the familiar 
mode functions of the Minkowski space time. The wavelength of a given mode becomes equal to the 
horizon size (or "crosses" the horizon) when kr] = 1. Inflation proceeds with r/ ^ 0, so the modes with 
progressively larger k cross the horizon. After horizon crossing, when ki] <^ I, the asymptotics of mode 
functions are 

i Uk iH 

The field variance is given by 

/d?k 
J^nf Iv^fcl' • (101) 

and we find in the asymptotic (the careful reader will recognize that this is already regularized expression 
with zero-point fluctuations being subtracted) 

(2^ J T 

Defining the power spectrum of the field fluctuations as a power per decade, (99^) = f Pip{k) din k, we 
find ^ 

P^{k) = ■ (103) 



i^') = 7KZ^i -■ (102) 



Curvature perturbations. According to Eq. the three-dimensional curvature of space sections of 
constant time is inversely proportional to the scale factor squared, ^^'^R oc a^^. Therefore, the perturba- 
tion of spatial curvature is proportional to 6a /a, and this ratio can be evaluated as 

C^- = F5t = i/^. (104) 

a If 

This allows to relate the power spectrum of curvature perturbations to the power spectrum of field fluc- 
tuations 

P^{k) = — P^{k) , (105) 
and we find for the power spectrum of curvature perturbations 

^C(^) = i ^ • (106) 

This very important relation describes inflationary creation of primordial perturbations, and can be con- 
fronted with observations. The usefulness of curvature perturbations for this procedure can be appreci- 
ated in the following way: 

1. Consider the perturbed metric, Eq. d46b . The product a(l — for the long-wavelength per- 
turbations can be viewed as a perturbed scale factor, i.e. 5a/ a = Comparing this relation with 
Eq. (11041) and Eq. d50b . we find for the temperature fluctuations which are of the superhorizon size at the 
surface of last scattering 

6T 2 

^ = 3 Cfe • (107) 

2. On superhorizon scales the curvature perturbations do not evolve.^ This fact allows to relate 
directly the observed power spectrum of temperature fluctuations to the power spectrum of curvature 
fluctuations generated during inflation. 

Tensor perturbations. Mode functions of gravity waves (after rescaling by Mp\/\/327r) obey the same 
equation as mode functions of massless minimally coupled scalar |92|. Using the result Eq. (I103t we 
immediately find L97J 

PH.) = 2 |f P.Ak) = 1^ 1^ , (108) 
where the factor of 2 accounts for two graviton polarizations. 

Slow-roll approximation. During inflation, the field ip rolls down the potential hill very slowly. A 
reasonable approximation to the dynamics is obtained by neglecting in the field equation if + 3Hip + 
V = 0. This procedure is called the slow-roll approximation 

<P»-^. (109) 

Field derivatives can also be neglected in the energy density of the inflaton field, p ^ V 

Svr 



t2 ^ 

This gives for curvature perturbations 



= . (110) 



^,^p(A;)l/2 = J^ = i^il. (Ill) 



should warn that this is quite a generic statement and does holds in situations usually considered. Thus, it is forgotten 
sometimes that this is not a universally true statement. To avoid possible confusion when encountering specific complicated 
models, the reader should keep this fact in mind. 



Normalizing to CMBR. As an example, let us consider the simplest model V = ^m'^ip'^. We have 



^ = ^, ^ndH=.^^. (112) 
y 2' V 3 Mpi ^ ' 



This gives for the curvature fluctuations 



IQiT niip 



Using the relation between curvature and temperature fluctuations, Eq. (11071) . and normalizing 5T/T to 
the measured value at largest /, which is 6T/T 10^^ (see Fig.|51l we find the restriction on the value 
of the inflaton mass in this model: 

^^^^..1013 GeV. (114) 
T 30 

Here I have used the fact that in this model the observable scales cross the horizon when (p « Mpi. 



Slow-roll parameters. The number of e-foldings (a = e^* = e^) of inflationary expansion from the 
time when (p = ipi to the end can be found as 

N{ifi) = / H{t)dt = -d^ = —^ —dip. (115) 

In particular, in the model Eq. (l94l we find that the largest observable scale had crossed the horizon {N ~ 
65) when ipi 3.5Mpi. All cosmological scales which fit within the observable universe encompass a 
small interval within Mpi < ip < ipi. And inflaton potential should be sufficiently flat over this 
range of Acp for the inflation to proceed. This means that observables essentially depend on the first 
few derivatives of V (in addition the the potential V{cl)o) itself). From the first two derivatives one can 
construct the following dimensionless combinations 



6=-^(-) , (116) 



which are often called the slow-roll parameters. 

The power spectra of curvature, Eq. ( I105t . and of tensor perturbations, Eq. (11081 . in slow-roll 
parameters can be rewritten as 

W = ^|^. ^-W = ^|^. (1.8) 



Comparing these two expressions we find 

Prik) 



16e . (119) 



Primordial spectrum. In general, the spectra can be approximated as power law functions in k: 

/ u \ "s-i 

P^ik) = P^iko) (-] , (120) 
Prik) = Priko) i^j • (121) 
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Scalar spectral index n^^ 

Fig. 20: Left panel: 95% constraints on inflationary models in the (n^, r) plane. From Ref. 1101 . Right panel: Forecast for 
the planned ESA Planck mission. 

To the first approximation, H in Eq. (II 18t is constant. Tlierefore, in this approximation, power spectra 
do not depend on k and ns = 1, ut = 0. This case is called the Harrison-Zel'dovich spectrum l98ll99l 
of primordial perturbations. However, in reality, H is changing, and in Eq dllSt for every k one should 
take the value of H at the moment when the relevant mode crosses horizon. In slow roll parameters one 
then finds (see e.g. Ref. LIOOJ for the nice overview) 

ns = l + 2r]- 6e, nr = -2e . (122) 

We can re-write Eq. (I123t as a relation between the slope of tensor perturbations and the ratio of power 
in tensor to curvature modes 

This is called the consistency relation to which (simple) inflationary models should obey. 

Different models of inflation have different values of slow-roll parameters ry and e, and therefore 
can be represented in the (77,6) parameter plane. Using the relations Eq. (11221) we see that this plane can 
be mapped into {ns, tit), or using also Eq. (I123t into the {ns, r) parameter plane, where r is the ratio 
of power in tensor to scalar (curvature) perturbations. In this way, different inflationary models can be 
linked to observations and constraints can be obtained. 

The most recent constraints in the {ns,r) plane, utilizing WMAP and SDSS data, are presented 
in Fig. |20l the left panel. The shaded dark red region is ruled out by WMAP alone. The shaded light 
red region is ruled out when adding SDSS information. The two dotted lines delimit the three classes of 
inflationary models known as the small-field, large-field and hybrid models. Some single-field models 
of inflation make highly specific predictions in this plane, as indicated. From top to bottom, the figure 
shows the predictions for V{(j)) oc cf)^ (line segment; ruled out by CMB alone), ¥{(()) cx 0^ (star; on 
verge of exclusion) and V{(j)) oc (jp' (line segment; the inflation model Eq. ( |94T ); still allowed). 

Testing inflation. All predictions of Inflationary cosmology, which could have being tested so far, 
have being confirmed. In particular, the Universe is spatially flat (within experimental errors), i7 = 
1.02±0.02, see Table I. The primordial perturbations are of superhorizon size and adiabatic. The spectral 
index is close to the Harrison-Zeldovich case, see Fig. |20| the left panel. Crucial test of inflationary 



models would be detection of gravity waves and verification of the consistency relation. This signatures 
of typical inflationary models are within reach of future CMBR experiments, see Fig.|20l the right panel. 

8. Ultra-High Energy Cosmic Rays 

In early years, cosmic ray studies were ahead of accelerator research, starting from the discovery of 
positrons, through muons, to that of pions and strange particles. Today we are facing the situation that 
the puzzling saga of cosmic rays of the highest energies may again unfold in the discovery of new physics, 
now beyond the Standard Model; or it may bring to life an "extreme" astrophysics. 

Immediately after the discovery of the relict Cosmic Microwave Background Radiation (CMBR), 
Greisen, Zatsepin and Kuzmin 11011 11021 have realized that the highest energy protons should catas- 
trophically loose energy in photo-production of pions on this universal background. This process limits 
the distance to the farthest sources of observed rays to be roughly 100 Mpc and should lead to the cut-off 
in the energy spectrum. However, the number of events with energies beyond the expected cut-off as 
measured by different installations is growing with time 1 103^ fT04l IT05l ITMl fT07 108 109 1, while no 
nearby sources where identified. The findings of Greisen, Zatsepin and Kuzmin (GZK) are based on 
solid fundamental physics which involve precisely measured cross-sections in a GeV energy range (in 
the center of mass reference frame). Therefore, if the data are correct - and it is believed they are basi- 
cally correct^ - one should either invoke new physics, or accept unconventional and uncomfortable very 
"extreme" astrophysics. This is the reason for the excitement and growing interest in ultra-high energy 
cosmic ray research; for recent reviews see LllO.. 111. .112111131 . 

Methods of detection. At energies below 10^^ eV, the flux of cosmic rays is sufficiently high that 
direct measurements using high altitude balloons or satellite experiments are possible. Above 10^^ eV, 
the flux is only one particle per m^ per year, which excludes direct observations on the orbit. At 10^'^ 
eV the number is down to one particle per square kilometer per century. Here the problem for direct 
measurements would be not only a vanishingly small flux, but the enormously high energy itself. (Re- 
member that calorimeters at modem colliders weigh hundreds of thousands of tonnes.) Fortunately, the 
major part of our UHECR detectors is already built for us by Nature and is put, rotating, into orbit: the 
Earth's atmosphere makes a perfect calorimeter. The atmosphere is just thick enough so that showers of 
secondary particles produced by incoming cosmic rays of the highest energies, in collisions with nuclei 
of air, reach their maximum intensity just above the Earth's surface. Particles in a shower propagate with 
the velocity of light, forming a thin disk perpendicular to the direction of the incident particle. At 10^^ 
eV the footprint of the shower on the ground is several kilometers across. 

The shower can be registered either by placing an array of particle detectors on the earth's surface, 
or by measuring the Cherenkov light produced by particles in the atmosphere, or by tracking the fluo- 
rescence light emitted when shower particles excite nitrogen molecules of the air. Particle detectors in a 
ground array can be spaced hundreds of meters apart and are operational around the clock. Fluorescence 
light telescopes see the cosmic ray track just like a fly's eye would see the meteorite, but only moving 
with the speed of light. These detectors are operational only on clear moonless nights, but are able to 
measure the longitudinal shower profile and its maximum depth directly. 

With either technique, the energy and incident direction of primary particle can be measured 
shower by shower. Chemical composition also can be inferred, but only in a statistical sense, after 
averaging over many showers. 

1 . Arrival direction. The timing of a signal in different detectors is used to determine the direc- 
tion of a shower (ground array technique). Direction is measured with an accuracy of about 2°. The 
measurement is straightforward and does not involve any uncertainties. Inferred information is reliable. 

'while recently a disagreement in measured fluxes has emerged, there is no reason for doubt in the reality of super-GZK 
events. 



Fluorescence light telescopes observe the whole shower track, and in stereo mode the precision of angle 
determination is 0.5°. 

2. Energy. Energy estimate, on the other hand, is not that straightforward. In fluorescent light 
detectors, the energy of primary particle is derived from the observed light intensity, therefore incorrect 
modeling and/or insufficiently frequent monitoring of atmosphere transparency can be a source of errors. 
For the ground array detectors, the energy estimate relies on a Monte-Carlo model of shower development 
and is related to the shower density profile. Nevertheless, the currently favored model, QGSJET II114I . 
describes data well from TeV up to highest energies and it is believed that the overall error (statistical 
plus systematic) in energy determination does not exceed 30 %. 

The best would be to employ both the ground array and fluorescent light techniques simultane- 
ously. This should reduce systematic errors, and this is the design of the forthcoming Pierre Auger 
project [?]. 

3. Chemical composition. Chemical composition can be inferred from the details of shower 
development. For example, showers initiated by heavy nuclei start earlier in the atmosphere and have less 
fluctuations compared to proton showers. Fluorescence detectors observe shower development directly. 
Using ground array detectors, the shower depth can be extracted by measuring e.g. the ratio of electrons 
to muons. At lowest energies, the chemical composition of cosmic rays reflects primary and secondary 
element abundances; for a recent review see 1 1 15 1; at highest energies, E > Ax 10^^ eV, the conclusion 
is that less than 50% of primary cosmic rays can be photonic at 95% confidence level L116J . 



8.1 Propagation of the ultra-high energy cosmic rays 

In this subsection, we consider the influence of different cosmological backgrounds on the propagation 
of highest-energy cosmic rays. 



Magnetic fields. Magnetic fields play an important role in the processes of cosmic ray acceleration and 
propagation, their trajectories being bent by the action of the Lorentz force 

dv Ze , ^, 

- = -[,TxB]. (124) 

For a qualitative discussion, it is often sufficient to compare a gyro-radius of the trajectory of a relativistic 
particle 

= ^ (125) 

to other relevant length scales of the problem. E.g., a magnetic "trap" can not confine a cosmic ray 
if the gyro-radius exceeds the trap size. The deflection angle A0, after traversing the distance L in a 
homogeneous magnetic field, is proportional to L/i?L- In a chaotic magnetic field, the deflection angle 
will grow as \/L. Let us estimate a typical deflection angle, L/Rg, of a charged UHE particle after 
traversing Galactic or extra-galactic magnetic fields. 

1 . In the Galactic magnetic field, for particles coming across the Galactic disc, we have 

^ . 2.5" _^ _^ , (126) 

Z E 3 nG 1.5 kpc 

where 3 fiG is magnitude of the regular magnetic field and 1.5 kpc is the width of the disc. We see that 
protons with E > 10^^ eV escape our Galactic disk easily. Protons of smaller energy are trapped and can 
escape the Galaxy only by diffusion and "leaking" from the boundary. Cosmic rays with E > 10^^ eV 
should be extra-galactic, if protons. Even if CRs would be all iron nuclei, at ^ > 2 • 10^^ eV cosmic rays 
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Fig. 21: Full sky map of deflection angles for UHECRs with energy 4 x 10^^ eV after traveling 100 Mpc in an extra-galactic 
magnetic field. The coordinate system is galactic, with the galactic anti-center in the middle of the map. Positions of identified 
clusters are marked using the locations of the corresponding halos in the simulation. The map is obtained in a magneto- 
hydrodynamical simulation of cosmic structure formation that correctly reproduces the positions and masses of known galaxy 
clusters in the Local Universe. From Ref. Ill 171 



should be extra-galactic, otherwise strong anisotropy in the direction of the Galactic disc would have 
been observed. 

2. Extra-galactic magnetic fields have not yet been measured, except for the central regions of 
galaxy clusters 1 1 18 1. However, there is an upper bound on their strength from the (absence of) Faraday 
rotation of polarized extra-galactic sources I119l ll20l. This translates to the upper bound on deflections 
in extra-galactic magnetic fields 

M „ 1020 eV B (LA)i/2 

— < 2.5° ^ — , (127) 

Z E 10-9 G 10 Mpc ' 

where A is the coherence length of an extra-galactic magnetic fields and is believed to satisfy A < 1 
Mpc. However, extragalactic fields are strongly inhomogeneous, with amplitude changing by orders of 
magnitude from clusters to filaments, and from filaments to voids. Deflections in some directions, which 
do not cross clusters and strong filaments, may be small, otherwise deflections can be very large. This 
situation cannot realistically be described by a mean field. 

Only recently have attempts been made to simulate UHECR propagation in a realistically struc- 



tured universe 1117111211 . Results of Ref. 11171 are shown in Fig. |2T1 Additional motivation for this 
simulation was to obtain, in constraint simulations of the Local Structure, a realistic map of expected de- 
flections, which would reflect the positions of known clusters. Such a map can be used in the analysis of 
cosmic ray arrival directions. Resulting deflections do not exceed the resolution of UHECR experiments 
over most of the sky. About an order of magnitude stronger deflections were obtained in Ref. |121 1. 
There are two possible reasons for disagreement. First, simulations of Ref. |121 1 were unconstrained 
and therefore do not reflect our concrete local neighborhood. Second, variable resolution of Ref. \ 
was better in cluster regions, which is a possible reason for the larger obtained dynamical range between 
fields in clusters and filaments. Since in both simulations the magnetic fields are normalized to typical 
values in the core of rich clusters, their values in the filaments will be very different, with larger fields 
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Fig. 22: Energy of protons as a function of the distance propagated in CMBR for three initial values of energy at the source, 
10^^ 10^^ and 10^° eV respectively. 



outside clusters in Ref. I121II . A work aimed to resolve these differences is in progress. For now, I adopt 
the results of Ref. 03 and conclude that arrival directions of UHECR should point back to the sources. 
Charged particle astronomy of UHECR is, in principle, possible. 

Interactions with cosmic radiation backgrounds. Ultra-high energy cosmic rays cannot propagate 
elastically in cosmic backgrounds. They have enough energy to produce secondary particles even in 
collisions with CMBR (important for proton primaries) or radio photons (important for UHECR photons) 
or infrared radiation (important for propagation of nuclei). Most important is the reaction of pion photo- 
production for protons (or neutrons) propagating in relic cosmic microwave background left over from 
the Hot Big Bang. For the threshold energy of this reaction we find, in the laboratory frame, 

(m„ + rn-Tr)^ — ml 

E,Up + l^N + 7r)= " . (128) 

2i/^(l — cos a) 

Note, that in the derivation of this relation, standard Lorentz kinematic and standard dispersion relation 
between particle energy and momentum, E'^ = k'^ + m?', are assumed. If any of these are violated, the 
threshold condition in a laboratory frame may look different. For the black body distribution of CMBR 
photons with temperature T = 2.7° K we find 

Eth^5x 10^^ eV . (129) 

This reaction has a large cross section, being the largest at the A resonance. At half-width of the reso- 
nance 

cr ~ 300 ^6 « 3 X lO-^^^cm^ . (130) 
Density of CMB radiation is n ~ ~ 400 cm . This corresponds to the mean free path: 

L„ = (o-n)"^ w 8 X 10^^ cm 2.7 Mpc (131) 

In each collision ^ 20% of energy is lost (which is the mass of a pion). Successive collisions rob protons 
of energy, which decreases exponentially. The distance over which energy decreases by one e-fold is 



called the attenuation length. At the threshold, Eq. (I129I) . the attenuation length is large, La ~ 10^ Mpc 
(being determined by other processes, see below.) With increasing energy, it rapidly decreases and at 
energies above the A resonance for typical CMBR photons, E ^ 5 ■ 10^^ eV, the attenuation length is 
La ~ 10 Mpc. It follows that the energy of protons drops below 10^*^ eV after it travels the distance of 
order 100 Mpc almost independently upon initial energy, see Fig. (l22l . We conclude that 

[ Al ] Protons detected with E > 10^*^ eV should have an origin within R < i?GZK = 100 Mpc. 

We will call the corresponding volume a GZK-sphere (or GZK-distance). 

The reaction p + 7 — > p + e+e~ is sub-dominant. While it has a smaller threshold (by a factor of 
2me/m7r ^ 10^^), it also has a smaller cross section. But, it becomes important at sub-GZK energies. 
Attenuation length for this reaction is 10^ Mpc - a noticeable and important effect. 

UHE photons loose energy in 7 + 7 — > e+e~. The threshold for the reaction with CMBR photons 
is smaller by a factor of 21711 /mT^nip ~ 10~^ compared to the GZK cutoff energy. The cross-section 
decreases fast with energy, a = a'^me/s^, where y/s is the CM energy and cjt ~ 10~^^cm^ is the 
Thomson cross-section. Therefore, attenuation length has a minimum at the pair production threshold. 
For CMBR photons, this occurs at w 2 • 10^^ eV and La ~ 10 kpc. The attenuation length increases 
with energy, reaching GZK distance roughly at £' lO^'' eV. Photons with even larger energies are able 
to penetrate even larger distances - and this is important for many models - but in this energy range, the 
main contribution comes from poorly known radio-background, which brings some uncertainty in the 
attenuation length of the highest energy photons. 

Heavy nuclei loose energy in photo-dissociation. Here, the main contribution comes from the 
infra-red background which is also poorly known. But again, at S 10^'^ eV the attenuation length is 
comparable to the GZK distance II122I . 

The cut-off. It is easy to understand why a sharp cut-off in the spectrum of protons should appear. This 
happens because the attenuation length decreases rapidly with increasing energy. Assume a power law 
injection spectrum for UHECR, Jin{E) oc E^"', and let n(r) be the density of sources. Fluxes from 
individual sources decrease as r^^, which is compensated by volume integration, r'^dr. Therefore, the 
total flux registered at energy E should grow in proportion to the upper limit of volume integration 



if the distribution of sources, n(r), does not depend on r. Here, R{E) corresponds to the attenuation 
length, i.e. the distance from which cosmic rays with energy E can reach us. The attenuation length of 
protons with E <b x 10^^ eV equals 10^ Mpc, while the attenuation length at £^ > 5 • 10^° eV is only 
10 Mpc. We conclude that 

[ A2 ] The drop influx by 2 orders of magnitude at GZK energy is expected if the distribution of sources 
is homogeneous . 

A word of caution is needed here. Transition in R{E) from sub-GZK to super-GZK regime is not 
instantaneous. Therefore, a particular value of the drop depends upon the shape of the injection spectrum, 
i.e. on the value of a, see e.g. Refs. fT23lfT24lfT25l . 

8.2 Generation of UHECR 

The origin of cosmic rays and/or their acceleration mechanisms have been a subject of debate for several 
decades. Particles can be accelerated either by astrophysical shock waves, or by electric fields. In either 
case, one can estimate the maximum energy; with optimistic assumptions, the final estimate is the same 
for both mechanisms. In practice, the maximum of energy is expected to be much lower. 




■R(E) 




(132) 



Fig. 23: Chandra telescope X-ray image of the nucleus, jet, and hot spot of Pictor A. From Ref. 11271 . 



1. Shock acceleration. Particles are accelerated stochastically while bouncing between shocks. 
Acceleration can continue only if particles remain confined within an accelerating region, in other words 
until gyro-radius, Eq. ( I125t . is smaller than the size of the region. This gives 

^max = ZeBL . (133) 

2. Acceleration by an electric field. The latter can be created by a rapidly rotating magnetized 
neutron star or black hole. If motion is relativistic, the generated electric field is of the same order as the 
magnetic field, and the difference in electric potentials is ~ (B x L). This, again, reproduces Eq. ( 11331 
for the maximum energy. 

Known astrophysical sources with {B x L) big enough to give E'max ~ 10^° eV are neutron stars, 
active galactic nuclei (AGN) and colliding galaxies. 

The central engine of an AGN is believed to be a super-massive black hole powered by matter 
accretion. AGNs have two jets (one of the jets may be invisible because of the Doppler effect) of rela- 
tivistic particles streaming in opposite directions. Interaction with the intergalactic medium terminates 
this motion and at the ends of jets the radio-lobes and hot-spots are formed, see Fig.|23l The acceleration 
of UHECR primaries may occur either near the black hole horizon (direct acceleration), or in hot spots 
of jets and radio-lobes (shock acceleration). The host of different AGNs is now classified in one unified 
scheme, for a review see LI 26.1 . Depending upon the angle between the jet axis and the line of sight we 
observe different types of AGN. A typical radio galaxy, showing two strong opposite jets, is observed 
at angles approximately perpendicular to the jet axis. An AGN is classified as a quasar if the angle is 
smaller than 30°. If we look along the jet axis (angle < 10°), i.e. directly into the barrel of the gun, we 
observe a blazar. 

It should be noted that not all radio-galaxies are the same. There are Fanaroff-Riley (FR) type 
I and type II galaxies (radio-loud AGNs), and Seyfert galaxies (radio-quiet AGNs). Both types of FR 
galaxies may be the sites of UHECR acceleration, but the hot spots in FR type II galaxies are considered 
to be most promising 1 128 1. It is believed that when observed along the jet axis, FR type II galaxies make 
a parent population of Highly Polarized Quasars (HPQ - subclass of blazars), while FR type I galaxies 
make a parent population of BL Lacertae objects (BL Lacs - another subclass of blazars). 

As an example, the X-ray image of the powerful FR-II radio galaxy Pictor A taken by Chandra 
observatory is shown in Fig. |23l Radio observations of jets have a long history. Recently, Chandra 
started to obtain high resolution X-ray maps of AGNs which, surprisingly, revealed very long coUimated 
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Fig. 24: Left panel: The energy spectrum of cosmic rays with a zenith angle up to 45° as measured by AGASA 11081 . 
Numbers near data points reflect the number of events in the respective energy bins. Right panel: A compilation of data from 
different experiments: Akeno (filled squares), AGASA (open squares), HiRes-I and HiRes-II (open and filled triangles), two 
Yakutsk sub-arrays of Cherenkov detectors (open and filled circles). From Ref. 11331 . 



X-ray jets. E.g. the distance from nucleus to the hot spot in Pictor A is at least 240 kpc. It is hard 
to explain such long jets as pure leptonic, and it is possible that the population of relativistic electrons 
responsible for the X-ray emission is a result of photo-pion production by UHE protons 1 127 1. 

Now, for any acceleration mechanism and independent of the actual acceleration site (i.e. be 
it either the AGN's black hole or any of the hot spots) the momentum of highest energy particles is 
expected to point in the direction of the jet 11291 . In other words, if AGNs are sources of UHECR, 
arrival directions of high energy cosmic rays may point back to a (subclass) of a blazar family. Such 
correlations were indeed observed II29III30lll3I[II32t with BL Lacertae objects. 

8.3 UHECR spectrum. 

The largest statistic of UHECR events has been accumulated for over 12 years of operation by the 
AGASA air shower array of surface particle detectors. The spectrum measured by the AGASA is shown 
in Fig. |24j left panel. The dotted curve represents the theoretical expectation for a homogeneous distri- 
bution of sources and proton primaries. This theoretical curve exhibits the GZK cut-off at E ^ 10^° eV. 
Remarkably, AGASA had detected 1 1 events with higher energy and the data show no hint for cut-off. 

It is hard to argue against the reality of these findings. AGASA exposure is under control, 
and the only issue is the energy determination. AGASA events have an accuracy of ±25% in event- 
reconstruction resolution and 18% in systematic errors around 10^*^ eV 1 108 1. Added in quadrature this 
gives RMS error of energy determination to be ±30%. More importantly, the probability of an upward 
fluctuation to 1.5 times the true energy is 2.8%. There are too many super-GZK events, and with this res- 
olution a spectrum with GZK cutoff cannot be transformed into an excess of post-GZK events assuming 
spillover 1 112|. 

Recently, the HiRes group has reported results obtained with a telescope which measures atmo- 
spheric fluorescence light. The energy spectrum is shown in Fig.|24]by triangles, the right panel. The 
spectrum is consistent with the GZK cut-off, and there are 2 events detected with E > 10^° eV. Sys- 
tematic error in energy measurement was estimated to be 21%, systematic error in the aperture is not yet 



clear. HiRes employs a relatively new technique, with the following issues usually cited for improve- 
ment: atmospheric attenuation corrections should be based on nightly measurements and not averages, 
better energy caUbration and aperture calculation are called for, see e.g. II112I . 

The Yakutsk group uses a hybrid detection method, combining a ground array of particle detec- 
tors with telescopes which are measuring Cherenkov light produced by a shower in the atmosphere. A 
recently reported 1 133J spectrum, derived from air Cherenkov light measurements, is shown in the right 
panel of Fig. |24]by circles. At the low energy end it agrees well with the Akeno spectrum, and at the 
high energy end it is consistent, within errors, with the AGASA spectrum. Yakutsk and AGASA disagree 
significantly with HiRes at ~ 10^^ eV were statistical errors are negligible, which points out to some 
systematics. AGASA and HiRes can be reconciled at E < 10^'^ by e.g. —15% and +15% respective 
shift of energy 11341 . The discrepancy between two experiments at £^ > 10^° after the shift is only 2a. 
However, even after these shifts there are still 9 events with E > 10^° in the combined data set. 

8.4 The Puzzle 

These measurements are regarded as a threefold puzzle, because contrary to assertions [Al], [A2] 

[ PI ] No candidate sources are found within the GZK distance in the directions of > 10^'^ eV 
events; 

[ P2 ] The AGASA spectrum does not exhibit the GZK cutoff. 
And finally we have the third puzzling question: 

[ P3 ] Which physical processes are capable of producing events with these enormous energies ? 

Conjectures. With the assumption that all three pieces of the puzzle, [PI] - [P3], are correct, the 
situation becomes desperate. There were no solutions suggested which would not invoke new physics 
beyond the standard model, or very speculative astrophysics. In addition, all models require fine tuning, 
and many do not really solve all three problems. It is not possible to consider here all the suggested 
solutions. Ignoring for now the problem [P2], the situation with [PI] and [P3] does not become easier, 
but we can now restrict ourselves to the discussion of astrophysical solutions only. 

Ignoring in addition [P3], the simplest suggestion is to assume very large extra-galactic magnetic 
fields, which would randomize UHECR trajectories. However, the results of Ref. Ill 171 do not support 
such a conjecture. As we have mentioned already, a consensus regarding EGMF fields has not yet been 
reached and in Ref. [121 1, much stronger extragalactic magnetic fields were advocated. Nevertheless, 
even in this case, the conclusion was that the condition of global isotropy of UHECR arrival directions 
requires the "local" value of magnetic field to be rather weak, B <, 0.1 ^uG, which, in turn, leads to 
a large number of UHECR sources in the GZK volume, n > 100; for similar limits on the number 
of sources see also Refs. [I35| \T3B ©81- These weak EGMFs of Refs. (Till IHl I rule out the 
possibility of a single powerful radio-galaxy, which happened to be nearby 1.1391 . or a gamma-ray burst 
scenario |[140|, as a potential sources of UHECR. 

Another suggested astrophysical scenario was a "dead quasars" model 11411 . This model assumes 
that quasars, powerful in the past, retain the possibility to accelerate to the highest energies even after 
the accretion of matter is exhausted and a quasar cease to be visible electromagnetically. However, the 
process of acceleration to the highest energies in compact sources is inevitably accompanied by a strong 
TeV emission 11421 . Recent results obtained by several TeV telescopes, in particular, non-observation 
of strong TeV sources, rule out the "dead quasar" model |143]. In addition, in Ref. 1 144] it was found 
that known quasar remnants are typically distributed too anisotropically to explain the isotropic ultra 
high energy cosmic ray flux except in the unrealistic case where extragalactic magnetic fields of 0.1 /xG 
extend over many Mpc. 

A possibility that ultra-high-energy events are due to iron nuclei accelerated from young, strongly 
magnetized neutron stars in relativistic MHD winds has also been suggested pi45^. However, with 



realistic parameters of Galactic magnetic field, even iron nuclei do not propagate diffusively within 
Galaxy, which disfavors this model 1. 146 J . 

Any observational clue? Many quite different models were suggested for the resolution of the GZK 
puzzle. The majority of suggested models, which we have no space to consider here, employ a new 
physics of one sort or another. (The reader may consult UHECR reviews cited at the beginning of this 
section, but I believe that a review which would cover all the suggested possibilities does not exist.) 
Instead, let us consider the question of whether or not there is already a clue in the data as to which 
model may be correct. Hints, and, in principle, critical signatures are given by: 

• Spectral shape. We do not yet have enough data at the highest energies to constrain models. 
Spectra below 10^'^ eV point to the AGN model of UHECR origin, with protons being primaries 
I1471ll48t . 

• Chemical composition. Again, not enough data at the highest energies. An analysis of Haverah 
Park data at lower energies shows that above 10^^ eV, less than 30% of the primary cosmic rays 
can be photons or iron nuclei at the 95% confidence level fl491. In other words, at least 70% 
should be protons. 

• Large-scale anisotropy. Gives strong signatures. Not observed, which is a hint by itself. Some im- 
plications we had considered already, and may add that the non-observation of anisotropy towards 
the Galactic center has a potential of ruling out the model of UHECR origin based on decays of 
super-heavy dark matter 1 150..15L.152..153J . 

• Small-scale clustering. This is an observed 11541 . reliable feature. (Errors in angle determination 
are definitely small.) It is already statistically significant. Therefore, below I shall concentrate on 
this signature. 

Small-scale clustering. It was observed by different installations that arrival directions of UHECR 
are too close to each other and this happens too often 1 155. . 154. 156.. 157 1. In particular, the AGASA 
collaboration has observed 6 doublets and 1 triplet of cosmic rays with E > 4 x 10^^ within 2.5°. 
The chance probability of observing just a triplet under an isotropic distribution is only 0.9% |154|. 
Statistical significance of these clusters in the AGASA data set was considered by several authors. In 
Ref. 1 157 1, an analysis based on the calculation of an angular autocorrelation function was employed, and 
the probability P = 3 x 10^^ of chance clustering was obtained. This includes the penalty for the choice 
of the energy cut, while the angular bin was chosen to be fixed at 2.5°, which is a value previously 
accepted by AGASA, being consistent with the angular resolution. In Ref. |158|, this analysis was 
repeated and confirmed. In addition, two more conservative estimates were done. In the first, the penalty 
factor for the adjustment of the angular bin size was added. This returns P = 3 x 10^'^. This is a valid 
procedure, but it misses prior information about the angular resolution of the installation. In the second 
estimate, the bin size was kept fixed, but the whole data set was divided in halves. The "original data 
set" 1 155 1 was used to justify the bin size of 2.5°, while clusters in it were removed for the subsequent 
evaluation of statistical significance. This procedure returned P = 8 x 10^^. Again, this is a valid 
approach too, and can be safely used with future large data sets. However, I'd like to stress that it is 
not an evaluation of the statistical significance of 6 doublets and 1 triplet. It is no wonder that a smaller 
data set has reduced statistical significance. Finally, in Ref. 1 159 1 it was found that the AGASA data set 
manifests a P ~ lO^'^ chance probability of clustering above background using independent statistics of 
(cos 6')[o°, 10°]- I fiiid this value, P 10^^, to be fair estimate of the current significance of clustering in 
the AGASA data. 

Note the following: if clusters are real and due to sources, the number of events in "physical" 
clusters should be Poisson distributed. Therefore, with the current low statistics, it is expected that 
roughly half of installations should observe significant clustering, while another half should not see it 



11571 . There is no clustering in the current HiRes data 11601 11591 . However, with the current statistics 
there is no contradiction yet I161[|159]| . 

The study of small-scale clustering is very important. If clusters are real and not a statistical 
fluctuation, then UHECR should point back to sources and UHECR astronomy is possible. Real sources 
should be behind the clusters and the correlation studies make sense. Pursuing this strategy, one should be 
restricted to astrophysical sources with physical conditions potentially suitable for particle acceleration 
to the highest energies. Active galactic nuclei (AGN) constitute a particularly attractive class of potential 
sources. As we have already discussed, if AGNs are sources, those which have jets directed along the 
line of sight, or blazars, should correlate with observed UHECR events. It is intriguing that statistically 
significant correlations of UHECR with BL Lacertae objects were found L129J . 



9. Conclusions 

Cosmology and astrophysics give us firm evidence that the standard model of particle physics is limited. 
The standard model fails to explain baryogenesis, does not contain non-baryonic dark matter and has 
no room for massive neutrino. We now know that dark energy also exists, but we do not know why it 
exists. There seems to be too many coincidences between numerical values of cosmological parameters 
which describe the matter and energy budget. Contributions of baryonic matter, non-baryonic matter and 
dark energy are almost equal at the present epoch, while they have seemingly unrelated origin and could 
differ by many orders of magnitude. Cosmology just became a precision science and is already full of 
surprises; we can expect even more exciting discoveries in the near future. 
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